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Abstract. Given a sequence of complex square matrices, an, consider 
the sequence of their partial products, defined by p„ = p„-ia„. What 
can be said about the asymptotics as n ^ cx) of the sequence f{pn), 
where / is a continuous function? A special case of our most general 
result addresses this question under the assumption that the matrices a„ 
are an h perturbation of a sequence of matrices with bounded partial 
products. We apply our theory to investigate the asymptotics of the 
approximants of continued fractions. In particular, when a continued 
fraction is Zi limit 1-periodic of elliptic or loxodromic type, we show 
that its sequence of approximants tends to a circle in C, or to a finite 
set of points lying on a circle. Our main theorem on such continued 
fractions unifies the treatment of the loxodromic and elliptic cases, which 
are convergent and divergent, respectively. When an approximating 
sequence tends to a circle, we obtain statistical information about the 
limiting distribution of the approximants. When the circle is the real 
line, the points are shown to have a Cauchy distribution with parameters 
given in terms of modifications of the original continued fraction. As an 
example of the general theory, a detailed study of a g-continued fraction 
in five complex variables is provided. The most general theorem in the 
paper holds in the context of Banach algebras. The theory is also applied 
to (r, s)-matrix continued fractions and recurrence sequences of Poincare 
type and compared with closely related literature. 



1. Introduction 
Consider the following recurrence: 

3 1 
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Taking l/oo to be and vice versa, then regardless of the initial (real) value 
of this sequence, it is an fact that the sequence is dense in R. Here is an 
illuminating proof. 

Fix xq and note that x„ = F^'^^xq), where F is the linear fractional 
transformation F{x) = (|x — and as usual F^") denotes the n-th 

composition of F with itself. 

Next recall the correspondence between 2x2 complex matrices and linear 
fractional transformations: let a given 2x2 matrix act on the complex 
variable z by 

a b\ , ^ az + h 
[z) := 



c dj ^ ' ' cz + d' 

so that matrix multiplication correspond to composition of linear fractional 
transformations . 

A simple calculation shows that the matrix for F is diagonalizable with 
diagonal elements a = 3/4 + i\/7/4 and (3 = 3/4 — is/?/ A. It follows that 
for some A G GL2(C), 



where X = a/(3 and r, s, t, and u are some complex constants. Notice that A 
is a number on the unit circle and is not a root of unity, so that A" is dense 
on the circle. The conclusion follows by noting that the linear fractional 
transformation 

rz + s 

z 1-^ 

tz + u 

is non-constant and thus a homeomorphism from C to C, and must take the 
unit circle to R U {oo}, since the values of the sequence Xn are real. 

Letting xo = 3/2 and iterating the definition of Xn gives that Xn is the 
n'th approximant of the continued fraction: 

(1.1) 3/2 , 



3/2 



3/2 



1 

3/2 



and thus one can write down the equality 
(1.2) R = 3/2 



3/2 



3/2 - 



1 

3/2 
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This is true so long as one interprets the "value" of the continued fraction 
as the set of limits of subsequences of its sequence of approximants. In this 
paper we generalize such equalities. 

From here on we employ space saving notations for continued fractions. 
In particular, an arbitrary finite continued fraction 

ai 

bo + 



bi + 



02 



b2 + 



as 



an 

bn-1 + -r- 
bn 



will be denoted by the expressions: bo + K^^i{ai\bi), 

, ai 02 an 

and 

bo+K^. 

1=1 bi 

If as n ^ oo this sequence of finite continued fractions converges to a limit 
in C (respectively C), then the continued fraction bo + K'^i{ai\bi) is said 
to converge in C (respectively C). The limit, when it exists, is called the 
classical limit of the continued fraction. 

Another motivating example of our work is the following theorem, one of 
the oldest in the analytic theory of continued fractions: 

Theorem 1. (Stern-Stolz, [SHlllSlllZlJ Let the complex sequence {bn} satisfy 
^ < cxD. Then 

oo 1 

6o + K - 

1=1 bi 

diverges. In fact, for p = 0, 1, 

lim P2n+p = Ap^C, lim Q2n+p = Bp £ C, 

n— »oo n— >oo 

and 

AiBo - AqBi = 1. 

Here Pi and Qi are respectively the numerator and denominator polyno- 
mials in the sequence {bn} obtained by simplifying the rational function one 
obtains by terminating the continued fraction after i terms. The Stern-Stolz 
theorem shows that all continued fractions of the general form described in 
the theorem tend to two different limits, respectively Aq/Bq, and Ai/Bi. 
(These limits depend on the continued fraction.) Although sometimes limits 
for continued fractions are taken to be in C, throughout we assume the limits 
for continued fractions are in C. The motivation for this is that continued 
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fractions can be viewed as the composition of linear fractional transforma- 
tions and such functions have C as their natural domain and codomain. 

Before leaving the Stern-Stolz theorem, we remark that although the the- 
orem is sometimes termed a "divergence theorem" , this terminology is a bit 
misleading; the theorem actually shows that although the continued frac- 
tions of this form diverge, they do so by converging to two limits {Aq/Bq and 
Ai/Bi) according to the parity of their approximant's index. Thus although 
these continued fractions diverge, they diverge in such a way that they can 
still be used computationally. 

A special case of the Stern-Stolz theorem is a fact about the famous 
Rogers- Ramanuj an continued fraction: 



(13) i + S. t ^ t 

^+1+ 1 + 1 + 1 + •••■ 

The Stern-Stolz theorem gives that for > 1 the even and odd approx- 
imants of (jl.3p tend to two limiting functions. To see this, observe that 
by the standard equivalence transformation for continued fractions, (jl.3p 
has the same classical approximants (see the beginning of section 3 for this 
terminology) as 

1111 11 

1 + 



l/q + l/q + l/g2 + l/q^--- + 1/g" + l/q^ +■■■' 

The Stern-Stolz theorem, however, does not apply to the following con- 
tinued fraction also considered by Ramanujan: 

, , -1 -1 -1 

(1.4) 



l + q + l + q^ + l + q^+---' 

Recently in [2j Andrews, Berndt, et al. proved a claim made by Ramanu- 
jan in his lost notebook ([32], p. 45) about (|1.4|) . To describe Ramanujan's 
claim, we first need some notation. 

Throughout take g G C with \q\ < 1. The following standard notation for 
g-products is employed: 

n-1 

(a)o := (a;g)o := 1, (a)„ := (a; g)„ := (1 - a g^), if n > 1, 

fc=o 

and 



{a;q)oo ■=Ylil- aq''), kl < 1- 



A:=0 



Set LV = e^'^^l'^ . Ramanujan's claim was that, for |g| < 1, 
(1.5) 

/I 1 1 1 \ ,/r!-c^"+l\ {q'-,q% 

lim - = —uj" ' ' 



^^oo \1 - 1 + g - 1 + g2 - . . . - 1 + gn + ay V^-'^""V ('J;?' 
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where 



n := 



1 — auj'^ (w^g, q) 



oo 



1 — am {ujq, q) 



oo 



Ramanujan's notation is confusing, but what his claim means is that the 
hmit exists as n — > oo in each of the three congruence classes modulo 3, and 
that the limit is given by the expression on the right side of (jl.Sp . Also, the 
appearance of the variable a in this formula is a bit of a red herring; from 
elementary properties of continued fractions, one can derive the result for 
general a from information about the approximants in the a = case. 

The continued fraction (jl.ip . the Stern-Stolz theorem, and (jl.4p are, in 
fact, examples of the same phenomenon. We define this phenomenon and 
investigate its implications. 

Now (jl.ip is different from the other two examples in that it has subse- 
quences of approximants tending to uncountably many limits. In a previ- 
ous paper [10], the authors presented a theory for dealing with finite sets 
of limits. Here continued fractions having finite sets of limits are unified 
with examples such as p.ip . Indeed, all of the examples above, including 
(jl.ip . are special cases of a general result on continued fractions (Theorem 
[7] below). More generally our results hold for a large class of sequences in 
Banach algebras. To deal with all of these cases we employ the notion of 
the sequential closure of a sequence. 

Definition. The sequential closure of a sequence in a metric space is 
the set of limits of convergent subsequences J3 The sequential closure of a 
sequence {sn}n>i is denoted by 

In this paper we study sequential closures in the specific context of se- 
quences of the form 



where Di are elements in a unital Banach algebra and / is a function with 
values in a metric space, often compact. Usually in this paper Di is a 
sequence of complex matrices. 

When working with infinite continued fractions we frequently omit the 
symbol (c) for the sequential closure. As with (II. 2p . equalities involving infi- 
nite continued fractions are to be understood as being between the sequential 
closure of the continued fraction and the set on the other side of the equation. 



Thus, for example, the sequence {1, 1, 1, ... } has sequential closure {1} although the 
set of limit (accumulation) points of the set of values of the sequence is empty. Note that 
in a survey paper describing some of the research in this paper, the authors previously 
used the phrase "limit set", unaware of the use of this phrase in the theories of discrete 
groups and dynamical systems. Also, the notion used here should not be confused with 
the sequential closure of a set, which also occurs in topology; thus, one should not confuse 
the sequential closure of a sequence with the sequential closure of its image. We thank 
Peter Loeb for the suggestion of the phrase "sequential closure" in the present context. 
More motivation is also given in the remark following Theorem 1101 
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This paper is divided into sections, which are now surveyed. The main 
results of section 2 are Theorems O [3l and |4] which are the most general 
results of the paper. Theorem [2] provides the underlying idea and is stated 
at the generality of Banach algebras. In section 2 we also discuss recent 
results of Beardon [5] , which apply hyperbolic geometry to study the con- 
vergence of compositions of Mobius maps. This approach yields short and 
geometrically motivated proofs of many results in the analytic convergence 
theory of continued fractions which generalize to compositions of higher di- 
mensional Mobius maps. Some of the results of [5] are related to ours in as 
much as they deal with generalizing the Stern-Stolz theorem. Indeed, one 
of the conclusions of Theorem [2] is similar to Theorem 4.2 from [5]. Overall, 
the main difference is that the theorems of [5], which generalize the Stern- 
Stolz theorem, do not generalize the particular conclusion of the Stern-Stolz 
theorem that the continued fraction's even and odd approximants tend to 
two different limits, and instead focus on the fact of divergence. The ap- 
proach of this present paper is to generalize the convergence of the even and 
odd subsequences in the Stern-Stolz theorem. Another difference is that the 
results here also give information about the sequences Pn and Q„ in the 
Stern-Stolz theorem. Recovering such information from the Mobius maps 
approach seems problematic. 

A special case of Theorem[2]is Theorem^ which is used to prove Theorem 
[71 which in turn gives detailed information about the sequential closures of 
continued fractions. This result is studied in detail in sections 3 and 4. 
Sections 5 and 6 use Theorem [H to study (r, s)-matrix continued fractions, 
and linear recurrences of Poincare type, respectively. 

Section 3 focuses on limit periodic continued fractions of elliptic and lox- 
odromic types. We discover a tapestry of results which weave together the 
sequential closure, modifications of the continued fraction, and the asymp- 
totics of the approximants of a large class of continued fractions (including 
many which represent naturally occurring special functions). Those of ellip- 
tic type do not converge, but we find that their sequential closures are well 
behaved, and that their approximants have nice asymptotics. Our results on 
the elliptic case add results to studies such as ^23j. The foundational result. 
Theorem [71 actually treats the loxodromic and elliptic cases simultaneously, 
evaluating the continued fraction when it can have one, finitely many, or 
uncountably many limits in a single formula. 

Section 3 also addresses the statistics of the sequential closure. In partic- 
ular, suppose a continued fraction has an infinite sequential closure. Then 
which points in the set have the "most" approximants tending to them, and 
which have the "fewest"? Thus for example, how can one describe the dis- 
tribution of the approximants of the continued fraction (|1.2p for M? These 
questions are answered by considering the geometry of associated linear 
fractional transformation. We show, for example, that when the sequential 
closure is M, the approximants follow a Cauchy distribution. Finally, sec- 
tion 3 considers the converse problem and gives an algorithm for finding a 
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subsequence of approximants tending to any given point in the sequential 
closure. 

Section 4 studies a non-trivial example of the theory. The section con- 
cerns a particular continued fraction with five parameters (and four degrees 
of freedom) which generalizes not only the identity (|1.2p above, but also the 
continued fraction ()1.4p . In fact, the main result, Theorem IIOI also gives 
the limit of the Roger s-Ramanuj an continued fraction (11. 3p for \q\ < 1 as a 
special case. Thus in one identity, we obtain the evaluation of a continued 
fraction when it has one, finitely many, or uncountably many limits. An ex- 
ample of this theorem is a perturbation of (jl.ip . Specifically, the sequential 
closure of the continued fraction 

1 1 1 1 

^^■^> g + 3/2- q^ + 3/2- q^ + 3/2- + 

where \q\ < 1 is complex, can be described exactly. In fact, this sequential 
closure is a circle on the Riemann sphere. (Thus as a consequence, when 
\q\ < 1 and q is real, (|1.6p always has sequential closure M.) Viewing this 
circle as a linear fractional transformation of the unit circle {2;GC:|2:| = 1}, 

az + b 
cz + a 

it transpires that the parameters a, 6, c, and d are special cases of the basic 
hypergeometric function 201, with complex base q. 

More generally, in (jl.6p if the numbers 1 and 3/2 are changed so that 
the limiting recurrence for the convergents of the continued fraction have 
distinct characteristic roots, there is a coherent formula, in terms of the 2<i>i 
basic hypergeometric function, for the sequential closure regardless of the 
nature of the roots. Indeed, cases in which the characteristic roots have 
a ratio that is a root of unity yield a continued fraction with a finite set 
of limits. In this circumstance, the approximants of the continued fraction 
converge to the set of limits when the approximants are taken in arithmetic 
progressions. 

Remarks and notation: 

(i) All sequential closure equalities in this paper arise from the situation 

lim d{sn,tn) = 

in some metric space {X, d) . Accordingly, it makes sense to define the equiva- 
lence relation x on sequences in X by {sn} s< {in} lim^^oo d{sn-, tn) = 
0. In this situation we refer to sequences {sn} and {tn} as being asymp- 
totic to each other. Abusing notation, we often write Sn x tn in place of 
{sn} ^ {tn}- More generally, we frequently write sequences without braces 
when it is clear from context that we are speaking of a sequence, and not the 
nth term. Identifying the constant sequence {L, L,L, . . .} with the value L, 
it is clear that the statements lim„^oo Xn = L and Xn ^ L are equivalent. 

(ii) It is a fact from general topology that given a compact topological 
space X and a Hausdorff space Y, then any continuous bijection g : X ^ Y 
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must be a homeomorphism and g and its inverse must both be uniformly 
continuous. Under these assumptions an immediate consequence for sequen- 
tial closures is: If {sn}n>i is a sequence with values in X, then (c){g{sn)) = 

9{©{Sn))- 

(iii) Another basic fact is that If {s„} and are two sequences in 

some metric space satisfying Sn x tn, then ©(«„) = ©(tn)- Additionally, 
if / is some uniformly continuous function, then the following sequence of 
implications holds: 

Sn-tn =^ fiSrr) X /(t„) =^ ©(/(s„)) = ©(/(tn)). 

2. ASYMPTOTICS AND SEQUENTIAL CLOSURES OF INFINITE PRODUCTS IN 

UNITAL BANACH ALGEBRAS 

The classic theorem on the convergence of infinite products of matrices 
seems to have been given first by Wedderburn [511 [52]. Wedderburn's the- 
orem is maybe not as well known as it deserves to be, perhaps because 
Wedderburn does not state it explicitly as a theorem, but rather gives in- 
equalities from which the convergence of infinite matrix products can be 
deduced under an li assumption. Wedderburn also provides the key in- 
equality for establishing the invertibility of the limit, but does not discuss 
this important application of his inequality. It is not hard to see that Wed- 
derburn's equations hold in any unital Banach algebra. This result will be 
employed in our first theorem, which gives asymptotics for oscillatory diver- 
gent infinite products in Banach algebras. This theorem is then applied to 
the Banach algebra Tld{C) of dxd matrices of complex numbers topologised 
using the /qo norm, denoted by || • ||. 

For any unital Banach algebra, let / denote the unit. When we use 
product notation for elements of a Banach algebra, or for matrices, the 
product is taken from left to right; thus 

n 

YlA^ := AiA2---An. 

1=1 

Theorems with products taken in the opposite order follow from the the- 
orems below by taking the products in the reverse order throughout the 
statements and proofs. 

Proposition 1 (Wedderburn). In a Banach algebra U, let I be the unit and 
Ai, i > 1, be a sequence. Then J2i>i W^iW < °o implies that nj>i(-^ + ^i) 
converges in U. Moreover, all the elements of the sequence I + Ai are 
invertible if and only if the limit nj>i(-^ + ^i) invertible. 

The following corollary provides a convenient estimate of the convergence 
rate of the product. 
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Corollary 1. Under the conditions of PropositionUi l^t L = Y\i^i{I + At) 
and Pm = YliLiil + Ai) . Then 

(2.1) \\L-Pm\\ <e^'>i 11^*11 -e^i<.<^\\M\ = o ( ^^H^'H ) ' 

\i>m / 

Proof. Put 

P^ = {I + A,){I + A2)---{I + Am), 

and 

Q^ = (l + ||^i||)(l + ||A2||)---(l + ||^m||). 
Expanding the product for Pm gives 

(2.2) Pm = / + ^ ^ An-^ + ^ ^ AnjAn2 

l<ni<m l<ni<n2<m 

~l~ ^ ^ AjiiAji2An^ 

I<ni<n2<n3<r?i 

+ ■ ■ ■ + ^ ^ Anj^An2 • • • ^nm-i ~l~ ^1^2 ' ' ' ^m- 

l<ni<ra2<---<'^m-i<'Ti 

Similarly, 

Qm = 1 + ^ ^ ll^nill ~l" ^ ^ ll^nill ||^n2ll 
l<ni<m l<ni<n2<m 

+ ^ ^ 1 1 An-^ 1 1 1 1 Ari2 1 1 1 1 + ^2 1 1 ' " ' 1 1 1 1 • 

I<ni<n2<n3<m 

Thus for m > k, 

(2.3) llPm-^'fcll < (Qm-Qfc), 

and 

(2.4) \\Pm - III < Q™ - 1 < e^">i - 1. 

From the standard condition for the convergence of infinite products of com- 
plex numbers, the convergence of X^„>il|^n|| implies the convergence of 
nj>i(l + ll^ill) s-iid thus that the sequence Q„ is Cauchy. Hence by (j2.3|) . 
Pm is also Cauchy, and so ni>i(-^ + exists. 

Recall that an element x in a Banach algebra is invertible if ||x — /|| < 1. 
For ni>i(-^ + ^0 to be invertible, it is obviously necessary that the elements 
of the sequence I + Ai be invertible. We show that this is sufficient. Since 
Yli>i \\Ai\\ < oo, there exists j G Z+ such that J2n>j \\^n\\ < log(2). Then 
(p:f| gives that 

||(/ + ^,-+i)---(/ + ^,-+™)-I|| <e^">^H^"l' -1. 

Letting m — > oo yields 

lim 11(1 + Aj+i) •••(/ + Aj+m) -I\\< e^">^ H^"" - 1 < e'°s(2) -1 = 1. 
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Hence liuim^ooil + ^j+i) •••(/ + Aj+rn) 
the left by the invertible elements I + Ai, 

Proof of Corollary. From Proposition [H 



is invertible. Multiplying this on 
1 ^ * ^ j gives the conclusion. □ 



\L Prr 



n(/+^,.)- n (^+^^ 



i>l 



l<i<m 



< 



n (^+^^ 



l<i<m 



< e 



o 



□ 



There have been a number of theorems recently on the convergence of 
matrix products, see [3l[5l[2l[8l[l2l[I7l[42l[48l[49]. Also closely related to 
Wedderburn's theorem are Theorems 3.7 and 3.8 of [5], originally given in 
[13j, which gives essentially the same result, restricted to SL2{C). There 
is also Theorem 6.1 of Borwein et al, [8], which is very similar in flavor 
to an earlier result given by Artzrouni [3]. In [9] Borwein et al. give a 
generalization of the theorem in [8] . 

Our focus here is on cases of divergence and our results concern finding 
asymptotics for the nth partial products. These in turn can be used to 
describe the sequential closures. 

Here we briefly discuss the important work of Kooman, |24l [25l I26j . 
Kooman makes an intensive study of the asymptotics of perturbed ma- 
trix products, recurrence sequences of Poincare type, and compositions of 
Mobius transformations. Kooman proves a number of different results on the 
asymptotics of recurrences, depending on various assumptions. In particu- 
lar, in some of his results, the li assumption is loosened, which we maintain 
throughout, and other theorems treat cases of equal roots of the character- 
istic equation for the recurrences. It should be noted in the present context, 
however, that Kooman's results on iterated Mobius transformations in [26] 
treat the case of the outer composition sequence, while continued fractions 
arise from the inner composition sequence, and thus his work on Mobius 
transformations doesn't seem to have immediate applications to continued 
fraction theory. (Recall that given a sequence of functions {/,} from a set 
to itself, the inner composition sequence is the sequence of functions {Fn}, 
where -Fn = /i ° /2 o ■ ■ ■ o /n and the outer composition sequence is the se- 
quence of functions where G„ = fn ° fn-i o • • • o fi.) Our results, 
which have a stricter assumption on the perturbation in the matrix product 
case, give detailed information about the relations between various associ- 
ated limits and the asymptotics for continued fractions and generalizations. 
Some of this information is probably lost under looser assumptions. But, we 
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think that a careful study of Kooman's results will yield new results on the 
convergence of continued fractions and generalizations. Note that Kooman 
did use his results on recurrences to solve a problem of Perron, see Chapter 
7 of [24]. In section 6 we compare the applications of the Theorem [2] below 
to Poincare type recurrences with the related results of Kooman and others. 

We set some further conventions and fix notation. Let G be a metric 
space, typically a subset of C^, where C is the Riemann sphere and g is 
some integer g > \. Here C is topologised with the chordal metric and 
the corresponding product metric is employed for C^. (This is defined by 
taking the maximum of the metrics of all the corresponding elements in two 
^-tuples.) Let / be a continuous function from a compact subset (to be 
specified) of a unital Banach algebra U, (usually 9Jl(i(C)) to G. Typically 
we do not distinguish different norms, the correct one being supplied from 
context. In a topological space the closure of a set S is denoted by S. 

Our first theorem is a perturbation result giving the asymptotics of di- 
vergent infinite products in a unital Banach algebra. Although we will only 
use a special case of this result, we believe the general result is of sufficient 
interest to warrant inclusion, especially since the proof of the general result 
requires no additional work. We denote elements of the Banach algebra by 
capital letters to suggest matrices, which is the case to which the result will 
be applied. 

Theorem 2. Suppose {Mi} and {Di} are sequences in a unital Banach 
algebra U such that the two sequences (for rj = 1 and rj = —1) 



(2.5) 



-1 



are bounded and {Di — Mi} £ /i(U), that is, 
(2.6) ^IIA -Mill < oo. 

i>l 

Lei e„, = X]i>„ IIA - Mi||. Then 

(n \ / " 

IlAj [jlMj 

exists and F is invertible if and only if Di is invertible for all i > 1. Also, 
(2.8) P-(U^)(U^^') =Oien)- 



\i=l 



',1=1 



As sequences 
(2.9) 



llDi>:FllMi 



i=l 



i=l 
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and moreover 



(2.10) 



0{en) 



YlDi-FllMi 

1=1 1=1 

More generally, let f be a continuous function from the domain 



Fl[Mr.n>h\u \J mA 



n>h 



for some integer h > 1, into a metric space G. Then the domain of f is 
compact in U and /(HILi -^«) ^ /(-^IlILi ^i)- Finally 

(2.11) © = © f[ M,^ , 

and 




(2.12) © / n^o =© ukii^' 



We do not assume compactness of G so it is possible that the equahties 
in the theorem are between empty sets. When G is compact these sets are 
clearly non-trivial. Note that the conditions of the theorem imply that all 
the elements Mj are invertible. 

An interesting special case of Theorem [2] is when the elements Mi are 
unitary matrices. In this paragraph the following matrix norm is used: 

1/2 

|M|I = I E I 

It is clear that ||M|| = ^/d when M is a d x d unitary matrix (for then 

||M|p = tr(MM ) = tr(/) = d), and thus the hypothesis on the sequence 
Mi is satisfied. More generally, one can assume that the sequence of matrices 
{Mi} are elements of some subgroup of GLd{C) that is conjugate to the 
unitary group. This case is important enough that we distinguish it in the 
following theorem. 

Theorem 3. Let {Mi} be a sequence of elements of a subgroup of GLd{C) 
that is conjugate to the unitary group. Then, if {Di} is a sequence GLa^C) 
and {Di — Mi} G li, all of the conclusions of Theorem\M hold. 

The special case of Theorem [2] that will be applied in the next section is 
U = 9Jlrf(C), Mi = M, where M be a diagonalizable complex matrix with 
eigenvalues on the unit circle. Since M is diagonalizable, put M = GEG~^. 
Then M^ = CE^C-^, and so it follows that ||M'=|| < ||C|| • \\C-'^\\ and 
IIM^^II is bounded for A; G Z. Thus the boundedness hypothesis is satisfied 
and Theorem [2] simplifies to the following. 
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Theorem 4. Under the above conditions, 

F = lim I TT A 1 M-"" 

\i=l / 

exists mSDT^(C) and det{F) ^ if and only if all the elements in the sequence 
{Di} are invertihle. Moreover, \\F — (HiLi A)-^~"il =0{£n), and 

(i) II nr=i A - FM^ W = 0(e„) . Thus ©(11^=1 A) = ©(FM"). 

(ii) Let f he a continuous function from the domain 

{FM^ : n > /i} U J 1 11^4' 

n>h U=l J 

for some integer h > 1, into a metric space G. Then the domain of f is com- 
pact and /(nr=i A) ^ /(FM"). Hence ©(/ (IlLi A)) = ©{f{FAn). 

Note that because M is diagonalizable, this theorem enables one to deter- 
mine the exact structure of the sequential closure using Pontryagin duality. 

A special case of Theorem H] is Proposition 1 of [10] in which the eigen- 
values are roots of unity. It considers asymptotics, but not the sequential 
closure. This special case is roughly equivalent to Theorem 1.1 of [43]. 

Proof of Theorem Observe that 



i-1 



i-1 



\i=l 



\i=l 



i=l 

' i-1 



3=0 



i-1 



i-1 



i-1 



i=l 



j=0 



j=0 



i-1 



i=l 

n 



j=Q 



i=l 

where / is the unit element and 

^i-l 



i-1 



Hence 



I A II < 



i-1 



• II A - Mi 



j=0 



j=0 



<C\\Di-Mi 



for some real absolute bound C. The second inequality followed from the 
boundedness assumption on the sequences (12. Sp . By (|2.6p it follows that 
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Yli>o W^iW < ^'^d SO by Proposition [H it follows that F exists and is 
invertible when the -D,- are invertible for i > 1. 



(2.13) 



lim 

n— >oo 



F 




Thus we have proved that 
1 

0. 



Again from the boundedness of the sequences in (j2.5p . 



(2.14) 

That is, 



lim 

n— >oo 



i=l 



i=l 



l[Di^Fl[Mi. 

i=l i=l 

Using this and the boundedness of the sequences in ()2.5p gives that the 
domain of / is compact. Thus / is not only continuous, but is uniformly 
continuous. This uniform continuity and (j2.14p give 



lim 



i=l 



\i=l 



0, 



and so 



The sequential closure equalities in the theorem follow from the third 
remark in the introduction and the error estimates follow from Corollary [1] 
and the boundedness assumption. □ 

We conclude this section by comparing these results to some of those 
from the recent paper [5], which mainly focuses on applying the hyperbolic 
geometry of Mobius maps to the convergence theory of continued fractions 
with complex elements. Consider the following two results from [5] that are 
closely related the results of this section: 

Theorem 5 (Theorem 4.2 of [5]). Suppose that G is a topological group 
whose topology is derived from a right-invariant metric ao, and that (G, do) 
is complete. Let /i, /2, • • • be any sequence of elements of G. Then, for each 
k, there is a neighborhood J\fk of fk such that if, for all j, gj £ Mj, then 
(51 ■ ■ ■ 5n)(/i ■ ■ ■ fn)~^ converges to some element h of G. 

The above theorem shares some of the structure of Theorem [2j In par- 
ticular it gives the existence of a limit similar to the limit F in Theorem [2j 
The hypotheses are quite different, however, and asymptotics are not given 
in Theorem [5j Also, sizes of the neighborhoods are not provided. 

For the following corollary, some definitions involving hyperbolic geometry 
are useful. A Mobius map acting on is a finite composition of maps each 
of which is an inversion or reflection in some A^— 1-dimensional hyperplane or 
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hypersphere in M^. The Mobius group acting on is the group generated 
by these inversions or reflections. The conformal Mobius group, denoted 
7W AT , is the subgroup of those maps that are orientation preserving, which 
means that they can be expressed as the composition of an even number of 
such inversions. See [5l[^. 

Corollary 2 (Corollary 4.3 of [5j). Let fi, f2, ■ ■ ■ be any sequence of Mobius 
maps. Then, for each k, there is a neighborhood Nt of fk such that if 
gj G J^j, j = 1,2, . . . , then there is some Mobius map h such that for all z, 
^■(Si ■ ■ ■ gn{z), hfi ■ ■ ■ fn)) ^ as n — > oo. In particular, for each point z, 
lim„ gi ■ ■ ■ gn{z) exists if and only if lim„ /i • • • fn{z) exists. 

The differences with our theorem are that the setting in Theorem [2] is 
more general and the sizes of the neighborhoods are not given in Corollary 
[2j However, in the case of complex Mobius maps, in [5] it is shown that the 
neighborhoods Nj can be taken to be the set of Mobius maps g that satisfy 

^^■^^^ " ^ ^^^TOF^WmP' 

Here the norms are of the matrix representations of the Mobius maps fi and 
9- 

Comparing this with Theorem [21 it can be seen that for the case of com- 
plex Mobius maps, unless enough of the norms are small, one expects 
our condition {Di — Mj} S li to be weaker in general, and thus our result to 
be stronger. Note that Theorem [2] also gives information about the sequen- 
tial closure as well as asymptotics with error terms. Information about the 
sequential closure is implicit, however, in Corollary 4.3 of [5j above. 

There is another theorem in [5] which is also related to Theorem [2j In 
fact, it is a generalization of the Stern-Stolz theorem presented in the intro- 
duction. Before stating the theorem, a couple of definitions concerning the 
hyperboloid model of hyperbolic space are required. 

For X and y in M^"*"^, let 

q{x, y) = xiyi + X2y2 H V xnVn - xn+iVn+i, 

and 

nN = {x e M^+^ : q{x,x) = 1, xn+i > 0}. 

TC]\f IS one branch of a hyperboloid of two sheets. It can be shown that 
TC]\f can be endowed with a hyperbolic metric and that the matrix group 
0~^{N + 1, 1), which preserve q as well as the condition xn+i > 0, act as 
isometries on this space. Let 5 be a Mobius map which acts on M^, and 
hence by the Poincare extension, on H^"'"^. Suppose then that g corresponds 
to the (A^ + 2) X (A^ + 2) matrix A which acts on Hn+2- In [S] the following 
beautiful generalization of the Stern-Stolz theorem is given: 

Theorem 6 ("The General Stern-Stolz Theorem" [5j). Suppose that gi, g2 
. . . are Mobius maps in M.n , and that gn is represented by the {N + 2) x 
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(A^ + 2) matrix An as above. If 



oo 



(2.16) 



T.VUn\\'-\W 



n=l 



converges, then the sequence gi - ■ ■ gn is strongly divergent. 

Consider the = case. Then this theorem should be compared with 
the case of Theorem [3] in which H is unitary, and the matrices M, represent 
Mobius maps. In Theorem [6l (j2.16p is exactly the condition required for 
Z^ra pO' to t)e bounded in M. (Here p is the hyperbolic metric on 
M, where M := {(xi,X2,X3) £ : X3 > 0}, and j = (0,0,1).) Next, 
(j2.16p is sufficient to guarantee that p{j,gi ■ • •(7n(j)) is finite, and thus the 
orbits of the product gi - ■ ■ gn never leave H. This later condition is what 
is meant by "strong divergence" . Now the Mobius maps that fix j are the 
unitary maps and g{i) = j if and only if H^lP = 2. The condition (j2.16p can 
thus be interpreted as saying that the elements gn approach some sequence 
of unitary elements sufficiently rapidly. This is roughly the same as the 
condition on the sequence {-Dj} in Theorem [3] when H is unitary. Of course 
the conclusion of the theorems go in different directions. 

In the next section we apply the d = 2 case of Theorem H] to get detailed 
information about the sequential closures of continued fractions. 



The signal result of this section is Theorem [7] which for the unifies the 
evaluation (in terms of sequential closures) of h limit 1-periodic continued 
fractions in loxodromic and elliptic cases. We also study the statistics of the 
classical approximants in cases of sequential closures of infinite cardinality. 

First recall the correspondence between 2x2 matrices and continued 
fractions. It is easily understood by remembering the correspondence be- 
tween compositions of linear fractional transformations and products of 2 x 2 
matrices, and then noting that the composition of linear fractional transfor- 
mations can be written as a continued fraction. For example, observe that 
for a general linear fractional transformation (avoiding cases such as c = 0): 



Thus, generically, any composition of a finite number of non-trivial linear 
fractional transformations can be written as a finite continued fraction. But 
to generate a continued fraction, one does not need to work with such general 
linear fractional transformations. In particular, working with transforma- 
tions of the form 



leads to the correspondence between matrices and continued fractions that 
will be used below: 



3. h LIMIT 1-PERIODIC CONTINUED FRACTIONS 




cz + d c - + z 
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(3.1) 
where 






Qn ° h + h + h + ' " + K' 



Here P„ and Qn are the numerator and denominator polynomials (called 
convergents) in the variables Oj and 6j obtained by simplifying the rational 
function that is the finite continued fraction. Their ratio, Pn/Qn, is called 
the nth classical approximant of the continued fraction. When working with 
continued fractions with arbitrary sequences of elements, it is assumed that 
at most one numerator element Oj vanishes. Also, we assume that for all n, 
{(in-,bn) / (0,0). Moreover, if for some positive integer A^, a^v = 0, then the 
continued fraction is finite, and thus converges. Notice that in this case, the 
sequences Pn and Qn still exist and are in general non-constant for values 
oi n > N . However, in this situation the ratio Pn/Qn will be fixed, when 
it is defined, and in the cases considered, it will be proved that the ratio is 
defined for n sufficiently large. 

From ()3.ip one reads off immediately the fundamental recurrences for the 
convergents P„ and Qn'- 



(3.2) 



Pn Pn-1 \ ( Pn-1 Pn~2 \ / &n 1 



^Qn Qn—1 / \Qn—l Qn—2 , 

/ \ ttn 



Taking the determinant on both sides of (|3.ip gives the determinant formula 
for the convergents of a continued fraction: 

(3.3) PnQn-i - Pn-iQn = (-l)"-'aia2 • • • a„. 

An infinite continued fraction 

(3 4) ^ On _ Oi 02 as 

n=l bn bl + b2 + b3 + ■ ■ ■ 

is said to converge in C (respectively in C) if 

ai a2 as an 
lim — — — 



n—foo 



bl + b2 + bs + ■■■ + br, 



exists in C (respectively in C). Let {un} be a sequence of complex numbers. 
If 

ai a2 as a„ 
hm — , — , — 



,.^oo bl + b2 + b3 + ■ ■ ■ + bn + UJn 

exist, then this limit is called the modified limit of K'^^i{an\bn) with respect 
to the sequence Detailed discussions of modified continued fractions, 

as well as further pointers to the literature, are given in [28j. Other treat- 
ments of the analytic theory of continued fractions can be found in Henrici 
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pj, Jones and Thron Perron [36j, and Wall [SD] . Note that by <^ 
and ([321), 

(3.5) 6o + 



6l + 62 + ^3 + ■ ■ ■ + &n + ti^n Qn+ i^nQn-l 

In the following theorem, the sequential closure of the sequence of approx- 
imants of a general class of continued fractions is computed. It transpires 
that the sequential closure is a circle (or a finite subset of a circle) on the 
Riemann sphere. The result was obtained by considering the problem of 
finding a uniform generalization of the examples in the introduction. 

Our theorem concerns the asymptotics of the continued fraction 

-aP + qi -af3 + q2 -af3 + qn 

(3.6) 



a + (3 + pi + a + P + P2 + ••• + a + P + pn 

where the sequences pn and qn approach in li and the constants a and /3 
are points in the complex plane. Specifically assume that 



(3.7) ^ \pn\ + \qn\ < 00. 

Let 



n=l 



Sn ■■= max ( ^ kil I , 

\i>n i>n / 



and put 

fn{w) 



-a(3 + qi -aP + q2 -aP + qn 



a + P + pi + a + P + P2 + ■■■ + a + P + pn + w' 



so that /„ := Pn/Qn = /n(0) is the sequence of classical approximants of 
the continued fraction (j3.6p . We follow the common convention in analysis 
of denoting the group of points on the unit circle by T, and we extend this 
notation in the following way. Define the unitary characteristic to be the 
map T;^ with domain C and codomain 2*^, given by T\ = ©(A") C C. From 
well-known facts it follows that, 



0, if|A|<l; 

T, if |A| = 1 and A is not a root of unity; 

|g27riA:/m ; Q < /c < m}, if A is a primitive mth root of unity; 

00, if |A| > 1. 



We identify the unitary characteristic with its set of values. 

Theorem 7. Throughout this theorem let {pn}n>i, {qn}n>i be complex se- 
quences satisfying (|3.7p ." a ^ P be complex numbers with a//3 G C. Consider 
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the following limits: 

(3.8) a = a{a,(3) = lim a""(P„ - /3P„_i), 

n— >oo 

b = b{a,l3) = - lim /3-"(P„ - aP„_i), 

71— >00 

c = c(a, /3) = lim a""((5„ - (3Qn-i), 

n—>oo 

d = d{a,(3) = - lim ^""■{Qn - aQn-i). 

n— >oo 

When \a\ = 7^ 0, all four limits exist, b{a,f3) = — a(/?, a), and d{a,(3) = 
—c{P,a); the limits for a and c exist if \a\ > while the limits for b and 
d exist if \a\ < \(3\ . 
We have 

(3 9) + gi -«/? + 92 -a/3 + qs ^ a^a/p + b 

a + l3 + pi + a + f3 + p2 + a + f5 + p2,+--- cT^jp + d' 

(Possibly non-existent limits are always annihilated by Ta/(}-) Assuming 
\a\ = |/3| ^ 0, and that qn 7^ afi, for n>l, 

az + b 



(3.10) X /i((q//3)"+^) where h{z) 



cz + d^ 



(3.11) det(/i) = ad - 6c = (/? - a) Jl ( 1 - ^ I . 

n=l ^ 

Moreover, when \a\ = |/?| 7^ z/ either |c| 7^ or /or < n < m , 
ca" + 7^ w/ien |c| = and a/ (3 is a root of unity, then as n ^ 00, 

(3.12) |/„-M(«//3)"+')|=0(e„). 

Finally, when \a\ = = 1, the following asymptotics for the convergents 

Pn and Qn hold as n ^ 00 

(3.13) 



aa" + 6/3" 



a — /3 



0(e„) and 



a — /? 



0(en) 



Note that the right hand side of (j3.9|) stands for the set {{az + b)/{cz + 
(i)|2; S Ta//3} and reduces to a point when a/(3 is or cxd. 

Definition. The cardinality of the sequential closure of a continued frac- 
tion is called the rank of the continued fraction. Clearly this definition 
makes sense for any sequence in a topological space. 

Thus for the complex continued fractions covered by Theorem [71 the rank 
belongs to the set Z"'"U{c}. (Here c denotes the cardinality of M.) In general, 
by Bernoulli's theorem on continued fractions, complex continued fractions 
can also have rank Hq. Notice that ()3.9p gives a unified evaluation of li limit 
1-periodic continued fractions in both the loxodromic and elliptic cases. 

Remark. The identity (13. 9p holds in the case where the continued frac- 
tion terminates, that is when qn = a(3, for some n, say n = N. 
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Theorem [7] is foundational for what follows. We give two corollaries be- 
fore the proofs (of Theorem [7] and its corollaries). Further results follow 
the proofs. The next corollary gives enough information to identify (up to 
sign) the specific coefficients in the linear fractional transformation h in the 
theorem in terms of modifications of the original continued fraction. The 
succeeding corollary makes that identification. 

Corollary 3. Under the conditions of the theorem the following identities 
involving modified versions of (13. 6p hold in C. When \a\ > 



-af] + Qn-i -af3 + Qn 
+ a + (3 + pn-i+ a+pn 



(3.14) h{oo) = -= lim fn{-f3) 

C n—>oo 

-aP + qi -aP + q2 
= lim 

n^oo a + P + pi + a + p + P2 + 



-af3 + Qn-i -a(3 + gn _ 
+ a + l3 + pn-i+ 13 + Pn 

and for k £ 7j when \a\ = |/?| ^ 0, we have 



-a/3 + Qn 
+ a + f3+Pn+ UJn-k ' 

where 

LJn = 1 — T £ ^) n G Z. 

It is also possible to derive convergent continued fractions which have the 
same limit as the modified continued fractions in Corollary [3l The key is to 
simply transform them via the Bauer-Muir transformation, see |28) . As this 
will be used in the sequel, it is presented here. 

Definition The Bauer-Muir transform of a continued fraction bo + 
K(an|^ra) with respect to the sequence {wn} from C is the continued frac- 
tion do + K(c„|(i„) whose canonical numerators Cn and denominators 
(convergents) are given by C_i = 1, = 0, Cn = An + An-iWn, 

Dn = Bn + Bn-iWn for n = 0,1,2,..., where {An} and {Bn} are the 
canonical numerator and denominator convergents of bo + K(a„|6„). 

Thus the Bauer-Muir transformation gives a continued fraction whose nth 
classical approximant is equal to the nth modified approximant of a given 
continued fraction. 



When \a\ < \P\, 

(3.15) h{0) = -= lim /„(-«) 

-aP + qi -al3 + q2 
= hm 

ji-^oo a + P + pi + a + P + p2 + 



(3.16) 



lim 



-aP -\- qi 



lim fn{iOn-k) 
n— >oo 

-aP + q2 



*(x> a + P + pi + a + P + P2 + ■■■ 
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Proposition 2. [HESIEI] The Bauer- Muir transform of bQ + 'K{an\bn) with 
respect to {wn} from C exists if and only if 

An = On - Wn-l{hn + Wn) / 

for n = 1,2,3, ... . If it exists, then it is given by 

, ^ Ai C2 C3 

Oi + tDi + (12 + "3 + ■ ■ ■ 

where Cn = a„_is„_i, (i„ = 6„ + - Wn-2'Sn-i, a^^c? s„ = A„+i/A„. 

Because of the generality of the continued fractions in Corollary [3l no 
substantial simplification occurs when the Bauer-Muir transformation is ap- 
plied, so we do not present the transformed versions of the continued frac- 
tions in the corollary. 

The following corollary gives (up to a factor of ±1) the numbers o, b, 
c, and d in terms of the (convergent) modified continued fractions given in 
Corollary [3l 

Corollary 4. When \a\ = \(3\ ^ the linear fractional transformation h{z) 
defined in Theorem has the following expression 

_ A{C -B)z + B{A-C) 
{C-B)z + A-C ' 

where A = h{oo), B = h{0), and C = h{l). Moreover, the constants a, b, c, 
and d in the theorem have the following formulas 

a = sA{C-B), b = sB{A-C), c = s{C-B), d = s{A-C), 

where 



\ {A-B){C-A){B-C) 



It is interesting to note that the sequence of modifications of (j3.6p oc- 
curring in (|3.16p converge exactly to the sequence /i((a//5)"''~^) which is 
asymptotic to the approximants of (|3.6p . 

Dividing through the numerator and denominator of the definition of ojn 
by gives that the sequence w„ occurring in (|3.16p is either a discrete 

or a dense set of points on the line 

-aT + f3 
T + 1 ' 

according to whether a//3 is a root of unity or not. Observe that —uJn+2 is 
the nth approximant of the continued fraction 

—a/3 —a(3 —a/3 
"^^^ a + /3 + a + /3 + ••• + a + /3' 

which, except for the initial a + P, is the non-perturbed version of the 
continued fraction under study. That the sequential closure of cOn lies on a 
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line follows from Theorem [8] below. Combining the continued fraction for 
LJn with (13.16P and Theorem [7] yields the intriguing asymptotic as A; ^ oo: 

-a/3 + qi -a(5 + q2 -afi + 

a + P + pi + a + P + P2 + ••• + a + P + pk 
-aP + qi -aP + q2 



lim 



a + P + pi + a + P + P2 + • 



(3.17) 



-aP- 



'Jn-l 



+ a + P+Pn-l + 



-aP + q„ 

Pn 



-aP 



-aP 



-aP 



a + P + a + p-\ + a + P 



n — k — I terms 



The continued fraction on the left hand side is divergent, while its trans- 
formed version on the right hand side asymptotically approaches the k-th 
approximant of the continued fraction on the left as A; — > co. The relation at 
(jS.lTp is valid when ^ \pi\ + \qi\ < oo and \a\ = \P\ / 0, and can be viewed 
as a continued fraction manifestation of Theorem HI 



a + P + pn 1 
-aP + qn 



-1 



a 
p 



M 



1 



a + P 1 
-ap 



-1 



Proof of Theorem Define 

(3.18) n 
For later use, note that 

(3.19) M = ' ^ 
that for n € Z, 

(3.20) M 

and that for n G 

(3.21) 



where, to save space later, we have put Qn = (q^~"/3^~")/(/3 — a). 

Let Pn and Qn denote the nth numerator and denominator convergents of 
the continued fraction (j3.6p . By the correspondence between matrices and 
continued fractions (13.111 




(3.22) 






( 1\ 


1=1 


I 1 oj 



no. 
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Now assume \a\ = \P\ = 1. Clearly 

\\Dn - M\\oo = max{|p„|, 

and thus 

"^\\Dn-M\\oo < OO. 



n>l 



It follows that the matrix M and the matrices -D„ satisfy the conditions of 
Theorem m Thus there exists F £ 97t2(C) defined by 

(3.23) 



F = lim 

n— >oo 



= lim 
(3.24) 
= lim 

n^oo 

(3.25) 








1 


:) 


Pn 


Pn-l 


Qn 


Qn-1 



Pn Pn-l 
Qn Qn—lj 



1 



-Q 
1 









1 1 



-1 




a(3 



\ 



-^^1,1-^2,2 — 2-^2,1 — lim {PnQn-1 — Pn~lQn) 



aj3 



l—n al—n 



P — a 



J 



Taking determinants in (j3.24p gives an expression for det{F): 

1 _ ^ 

n=l 



The last equality follows from the determinant formula for continued frac- 
tions (I33|) . 

Consider the non-terminating elliptic case, i.e. \a\ = = 1 and g„ 7^ a/5 
for n > 1. It follows that det(F) / 0. Let / : GL2(C) ^ C be given by 



/ 



u V 

W X 



u 
w 



Note that / is continuous, and thus using TheoremHJ is uniformly continuous 
on the compact set 



{FM" : n > 1} U y I 

n>l 



Pn Pn~l 
Qn Qn—1 
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Theorem H] and the matrix product representation of continued fractions 
then give that 

Hence using (|3.2U|) and the definition of /, 

.3 El^ ^ Fi,i(a"+i - /?"+^) + Fi,2(-a/? (a" - /?")) 

^' ' Qn ^i^2,i(a"+i -/3"+i) + F2,2(-a/?(a'^ -/?")) 



n+l 



(Fi,i-/?Fi,2)(^fj +(aFi,2-Fi, 

(F2,l-/?i^2,2)(f) +(aF2,2-i^2, 

/i((a//3)"+i), 



where 

(3.27) h{z) 



az + h 
cz + 



with o = Fi^i - /?Fi_2, ^ = a-F'1,2 - -P^i,!, c = F2,i - /?F2,2, = 0^2,2 - -^2,1, 
and Fi^j G C are the elements of F. The hmit expressions for o, 6, c, and 
d in the theorem follow by simplifying the constants in h defined here, and 
then using (j3.25p . Next notice that from (j3.25p . the elements of the matrix 
F are symmetric in a and (3. This along with the symmetry of P„ and Qn as 
well as the definitions of a, 6, c, and d implies that 6(a, (3) = —a{(3, a) and 
also that d{a,f3) = — c(/?, a). The limits ()3.8p are clearly invariant of the 
size of |a| (since P„ is a polynomial of degree n in a), so they all exist under 
just the assumption |a| = 7^ 0. Note that we can compactly express the 
definition of a, b, c, and d in the following matrix equation: 

c d J \F2,i F2,2 J \-P a 

The product formula for ad — be follow immediately by taking the determi- 
nant and using the expression for det(-F) above. Solving for F gives 



(3.28) F 



a 



b \ /a 1 \ 1 



c d J \/3 1 J a-p 



Now /i : C ^ C is a continuous bijection when det(/i) 7^ 0. Put A = a/ (3. 
From (|3.26p and the remarks in the introduction, 



©(/i(A"+^)) = M(c)(A"+^)) = /i(T,), 

and so (j3.9|) is proved. 

To prove the asymptotics for P„ and Qn employ Theorem H] (i) to obtain 

(3.29) )xFAr, 
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Substituting ([M]) and IKm into (fM]) yields 

Pn Pn—1 
,Qn Qn—1 



a b 

C d 

a b 

c d 



aa 



ri+l 



ca 



,n+l 



+ bl3 
+ df3' 



a 1 

/3 1 

,n+l 
n+l 



a 



ri+l 
-n+l 



* 



-a(3{a'^ 
1 

a — P 
1 

a — f3 



n+l 



a/3" - /3a" 



1 



(a-/3)2 



Thus the sequences Pn and Q„ have the claimed asymptotics by Theorem 
Put ^„ = aa" + ^/3" and Bn = cd"- + d/?" and observe that 



|/n-MA"+')| 



< 







< 


Pn Bji Aji Bfi 


Qn 


Bn 




QnBn 


1 


en + 


A 


n 


e^n ) 


Qn 


QnBn 



+ 



-^nBn 



Qn-^n 



QnBn 



and this error is 0{en) providing that Bn is bounded away from 0. (Recall 
that Qn X Bn/{a — /3).) It is easy to see that Bn is bounded away from 
under exactly the two conditions given in the theorem. The restriction 
|a| = |/?| = 1 can be loosened to \a\ = |/3| / by employing the equivalence 
transformation in which the numerator elements of the continued fraction 
are divided through by |ap, and the denominator elements by |a|. 

Now consider the non-terminating loxodromic case , i.e. |a| ^ \(3\ and 
Qn 7^ a/3, and assume that a/3 ^ 0. Since the continued fraction (j3.6p is 
of loxodromic type, it converges by Theorem 28, p.p. 151-152 of [28j. By 
Theorem 5.1 of [24] the recurrence for the conver gents, 



Yn = {a + P + Pn)Yn-l + (-a/3 + qn)Yn-2 

has a basis of solutions {u„(a, /3), f„(a, /3)} satisfying 



lim 

n— >oo 



lim 

n— >oo 



1. 



a" n^oo /3" 

Thus there exist r(a,/3) and s(a,/3) such that 

Pn = r{a, P)un{a, /3) + s(a, /3)i;„(a, /3). 

By (|3.30p when |a| > |/3|, lim n — *cxD a ^ Pn — ^(a, /3); similarly lim^ — >oo^ ^ 
Qn exists. Hence the limits for the constants a and c in (jS.Sp exist. Since 
|a| > |/3|, T^/p = oo. Thus the right hand side of ()3.9p reduces to a/c. Now, 



a Pn- PPn-1 

— = hm 

C n^oo - PQn-1 



hm /„(-/3). 
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It is well-known, see j28j p. 160, that liuin^oo fn{— P) = when \a\ > 
\P\ > 0, where /i is the classical limit of the continued fraction. The case 
1/3 1 > \a\ > follows by symmetry. 

Next consider the non-terminating loxodromic case where |a| > /3 = 0, 
with Qn ^ for n > 1. By taking an equivalence transformation, this is 
equivalent to the case a = 1, P = 0. Again by Theorem 5.1 of [24] the 
recurrence for the convergents, 

y„ = (1 +p„)Yn-l + g'n>"n-2 

has a solution {n„(a,/3)} satisfying lim„_^oo'Un = 1- By the Poincare- 
Perron theorem (see the beginning of section 6) there also exists a solution 
{vn{a,(3)} satisfying liuin-^ooVn+i/vn = 0. This implies that there exists a 
A; > such that t>„ 7^ for n > k. Select such a k. Then for n > k, 

hm — = lim • = U, 

n^oo Vk n^co Vk ffc+l Vn^l 

and thus it can be concluded that lim„^oo Vn = 0. Since P„ and Q„ are 
linear combinations of Un and f „ , it follows that the limits for a and c in ()3.8p 
reduce to a = lim^^oo Pn and c = lim„^oo Qn- It is known that when a = 1 
and (3 = 0, the continued fraction converges in C to lim„_>oo Pn/Qn = ci/c; 
see [28J, p.p. 15H52. When a = 1 and /? = the right hand side of (l3^ 
simplifies thus: 

olTo//? + ^ _ aco + b _ a 
c^a/p + d coo + d c ' 

which as has just been shown is the classical limit of the continued fraction. 

Now consider the terminating elliptic case, that is assume that |a| = 
\(5\ = 1, and that = a/5 for some > 1. Theorem |4] still gives that 
F exists and that the asymptotics (13.131) also hold. It follows that a, 6, c, 
and d exist, although ad — be = 0. It is easy to check that in this situation 
a = a*a-^PN-i, b = b*f3-^PN-i, c = a*a-^QN~i, and d = b*p-^QN^i, 
where the numbers a* and b* are defined by the limits 



hm a'HPk - PPk-i), 



(3.31) 6* = - hm P~\P^-aPU), 

k—>oo 

and is the A;th numerator convergent of the continued fraction 

(3.32) b^+K -«^ + «^+» 



i- 



--ia + (3 + pN+i 



(Note that we have use the easily proved identities PN+k = -ffc-P/v-i and 
QN+k = PkQN-i which follow from the assumption that = Oij3.) Because 
(j3.32p has no vanishing numerators, it follows that it is impossible for both 
a* and b* to vanish. Indeed, what we have already proved above applies to 
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the continued fraction ()3.32p (so by (13. lip a*d* — 6*c* / 0). Computing the 
right hand side of (|3.9p gives: 



aT 



"//9 



N, 



P, 



N~l 



a a 



P 



N-l 



Qn-iJ a*a-^T^/p + b*P 



-N 



Q 



N-l 



since at least one of the hmits a* or b* is non-zero. The condition that 
= = 1 can now be loosened to just \a\ = |/3| / by applying an 
equivalence transformation. 

Finally, consider the terminating loxodromic case, i.e. \a\ > and 
qn = a(3. As Tq/^ = cxd G C, we need to show that the limits a and 
c exist and that a/c = Pn-i/Qn-i- The proofs in the non-terminating 
loxodromic cases apply here as well and show that the limits a and c ex- 
ist. Moreover, a = a*a~^ P^-i, and c = a*a~^ Qn-i, where a* is as 
defined in the terminating elliptic case above. Again, the same argument 
as in the non-terminating loxodromic cases give that lim^^oo Oi~^Pk 
ists. Moreover, from the Poincare-Perron theorem on recurrences, see sec- 
tion 6, limfc^oo -Pfc+i/^fc exists, and is equal to either a or (5. Thus, since 
a, /? G C, there follows P^ ^ Q for k sufficiently large. Hence it follows that 
Pn-i/Qn-1 = a/c in this case. The case \a\ < \(3\ and qjy = a(3 follows by 
symmetry. □ 

Proof of Corollary O ()3.14p and (j3.15p follow immediately from the value 
of a modified continued fraction (j3.5p . with uj^ = —13 and w„ = —ol, respec- 
tively, and the limit expressions for a, 6, c, and d. 

Let / denote the function from the last proof. To get ()3.16p , observe that 



= / 



f{FM^) = f 



lim 

n— >oo 



Pn 
Qn 



Pn-l 

Qn-1 



lim 

n—foo 



P P 

Qn Qn—1 



1 I M~ 



/ 



/ 



lim 



Pn 
Qn 



Pn~l 
Qn-1 



a 



n—k—1 



\ 



n—k 



13' 



a" 



n—k—l 



a 



af3 



n-k+l 



a 



n-k+l 



al3 



lim 

n— >oo 



(a 



n—k—l 



'n—k—l 



)Pn - (a"-*= - /3"-^')^„-i 



(^^n-k-l 
Pn 



a 



n—k 



_ j^n-k 



P'^-nQn- 



lim 



a 



n—k—l 



n—k—l 



P. 



n-l 



Qn 



a 



a 



n—k—l 



l^n-k 
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= iim — — 

n^oo Qn + OJn-kQn-1 

where 



The result now follows from (13.51). □ 



Proof of Corollary^ The expression for h{z) follows immediately using al- 
gebra from (j3.14p . (jS.lSp . and (j3.16p with k = —1. The expressions for a, 
b, c, and d follow by using (jS.lip along with the fact that the coefficients in 
the two expressions for the linear fractional transformation must be equal 
up to a constant factor. □ 

Note that putting k = and A: = — 1 in (j3.26p gives the following identi- 
ties: 



(3.33) /i(A) 

Ml) ^ - 

One naturally wonders just how effectively the parameters a, b, c, and d 
in Theorem [7] can be computed. In the next section, a particular contin- 
ued fraction is considered which generalizes one of Ramanujan's, as well as 
(jl.ip . and these parameters explicitly are computed as well-behaved mero- 
morphic functions of the variables in the continued fraction. Thus, for the 
gr-continued fraction studied in the next section, the parameters can not only 
be computed, but also have nice formulas. 

An interesting special case of Theorem [7] occurs when a and /3 are distinct 
m-th roots of unity (m > 2). In this situation the continued fraction 

-a/3 -I- qi -a/3 + q2 -af5 + gs -a/3 + 
a + P + pi + a + P + P2 + a + /3 + p3 + a + /3 + p4 + ••• 

becomes limit periodic and the sequences of approximants in the m different 
arithmetic progressions modulo m converge. The corollary below, which is 
also proved in [10] , is an easy consequence of Theorem [71 Note that by 
Theorem[71 ()3.6p can also have a finite sequential closure in the more general 
case that a//3 is root of unity, a case not covered in the following corollary. 

Corollary 5. Let {pn}n>i, {Qn}n>i be complex sequences satisfying 

oo oo 
^|pn|<00, \qn\ < OO. 

n=l n=l 

Let a and (5 be distinct roots of unity and let m be the least positive integer 
such that a™ = /3"^ = 1 . Define 

-a(3 + qi -a/3 + q2 -a/3 + 



G 



a + (3 + pi + a + /3 + p2 + a + (3 + pz + ■■■ 
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Let {Pn/Qn}'^=i denote the sequence of approximants of G. If Qn 7^ Oi(3 for 
any n > 1, then G does not converge. However, the sequences of numerators 
and denominators in each of the m arithmetic progressions modulo m do 
converge. More precisely, there exist complex numbers AQ,...,Am~i and 
Bq, . . . , Bm-i such that, for < i < m, 



(3.34) 



lim Pr^^fc^ 
K— ►oo 



lim Qmk^ 



Bi. 



Extend the sequences {Ai} and {Bi} over all integers by making them peri- 
odic modulo m so that Iji3.34\ ) continues to hold. Then for integers i, 



(3.35) 

and 
(3.36) 

Moreover, 
(3.37) 



B, 



AiBj 



Ai - PAq 
a — f3 

(3Bo 



B 



a — P 



a' + 



a' + 



aAo - Ai 

a- 13 

oBq — Bi 
a — P 



/3^ 



AjBi 



-(a/3) 



a 



•I- J 



a — P 



n=l 



aP 



Put a := exp(27ria/m), P := exp(27ri6/m), < a < 6 < m, and r := 
m/gcd(6 — a,m). Then G has rank r and its sequential closure is the finite 
set in C given by {Aj/Bj : 1 < j < r}. Finally, for k >0 and 1 < j < r, 

Aj-\-kr Aj 



I^j+kr 



B, 



Proof. Let M be as in Theorem [71 It follows from (j3.19p that 



(3.38) 

and thus that 

M™ = 



1 
1 



a — P 

aP {a^ -P^) 
a — P 



a-' 



a — P 

-a^ P + aP^ 
a-p / 



1 
1 



1 < j < m. 



Let the matrix F be as in Theorem [71 From (j3.29p . 



(3.39) 



lim 

n— >oo 



Q, 



rnn+t 



Qmn+i—1. 



lim FM™"+* = FM' 

n— >oo 



This proves ()3.34p . 
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Now let Ai := limn^oo Pmn+i, and Bi := limn^oo Qmn+i- Notice by 
definition that the sequences {Ai} and {Bi} are periodic modulo m. It 



easily follows from (I3.39[) that 



Ai j4j_i\ / Aj Aj—i 



Bi Bi-iJ \ Bj Bj^i^ 



(|3.38p also gives that 

(3.40) A = A,°'""-r"-^,_.°^(°"'-^'") 

a — p a — jj 

and 

(3.41) = 



Thus 



a — (3 a — (3 

{Aj - A^i+j Bj) a/3 {a^-^ - (3'-^) 



a — P 

Equations ([3:35]) and (f336|) follow from (fOO]) and ^^JT\i by setting j = 1. 
()3.37p follows after applying the determinant formula 

mk+j 

AjBj^i - Aj^iBj = - lim TT (a/3 - g„) 



n=l 
oo 



("«'n(i-||) 

n=l ^ ^ 



Since X^^i \ Qj\ converges to a finite value, the infinite product on the right 
side converges. 

For the continued fraction to converge, AiBi-i — Ai-iBi = is required. 
However, (13.370 shows that this is not the case. □ 

3.1. Distribution of approximants. Let T' denote the image of T under 
h, that is, the sequential closure of the sequence {/«}. The asymptotic for 
fn given in Theorem [7| is 

(3.42) /„ X /i(A"+i), 

where h is the linear fractional transformation defined in the theorem and 
A = a/(3. 

Some observations can immediately be made. It is well known that when 
A is not a root of unity, A"""*"^ is uniformly distributed on T. However, the 
linear fractional transformation h stretches and compresses arcs of the circle 
T, so that the distribution of h{X'^~^^) in arcs of T' is no longer uniform. 
(Recall uniform distribution on a curve happens when as n — > oo each 
segment of the curve get's the proportion of the first n points equal to the 
ratio of the segment's length to the length of the whole curve.) Additionally, 
T' may not be compact in C. So we consider a probability measure on T' 
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giving the probability of an element /i(A"'^^) being contained in a subset 
of T'. This measure is easy to write down. Let S C T', then h~^{S) is 
a subset of the unit circle. Then since A"' is uniformly distributed on T, 
P{S) := n{h-^{S))/2Tr gives the probability that for any n, h{X^) G S. 
Here n denotes the Lebesgue measure on T. Note that P depends entirely 
on h, and thus only on the parameters a, b, c, and d. 

In general /„ ^ T', but because of (|3.42p . as n ^ oo, the terms of the 
sequence fn get closer and closer to the sequence h{X^~^^) which lies on T'. 
Thus we speak of P as the limiting probability measure for the sequence fn 
with respect to T'. When the sequential closure is M, we speak of the limiting 
probability density function. 

More specifically, (j3.42p implies that there is a one-to-one correspondence 
between the convergent subsequences of /i(A"^^) and those of fn such that 
the corresponding subsequences tend to the same limit. As /i is a homeo- 
morphism and A" is uniformly distributed on T, it follows that the proba- 
bility of an element of ©{fn) being contained in a subset S of T' is exactly 
P{S)=fi{h-\S))/27r. 

Fortunately, this distribution is completely controlled by the known pa- 
rameters a, b, c, and d. The following theorem gives the points on the 
sequential closures whose neighborhood arcs have the greatest and least 
concentrations of approximants. 

Theorem 8. When a/ (3 = X €z T is not a root of unity and cd ^ 0, the 
points on 

aT + b 
cT + d 

with the highest and lowest concentrations of approximants are 

Oil a, I 6, ,1 

-|c| + -|d| — |c| + -|d| 

^ , and \ , , 

\c\ + \d\ -\c\ + \d\ 

respectively. If either c = or d = 0, then all points on the sequential 
closure have the same concentration. The radius of the sequential closure 
circle in C is 

a — (3 



n 1 



|c|2 - |(i|2 11 V aP 

n=l 

and its center is the complex point 

\hil)\'{hi-l) - h{i)) + \h{~l)\'m) -HI)) + \h(i)\Hh{l) - hj-l)) 
hil)ihii) - h{-l)) + h{-l)ihil) - h{i)) + h{i){h{-l) - h{l)) 

The sequential closure is a line in C if and only if \c\ = \d\, and in this 
case the point of least concentration is oo. 

Proof Let g{e) = h{e'^). Thus g{e) parametrizes T' for 9 G [0,2tt] and e'^ 
moves with a uniform speed around T as moves uniformly from to 27r. 
Then g{6) moves around T' at different speeds depending on how the length 
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g{9) change with 9. Accordingly, we wish to compute the rate of change of 
the length of g{9) with respect to 9. We then wish to know when this value 
is minimum and maximum. To this end put 

m := f \g'{t)\dt. 







Accordingly, I' {9) = \g'{9)\. An easy computation gives 

\ad — bc\ 



no) 

l"{9) = 



|c|2 + |(i|2 + cde*^ + cde-'' 

and thus _ 

\ad - bc\{cde-'^ - cde'^) 



(|c|2 + |d|2 + cde*^ + cde-^^)2 

Clearly l"{9) = if and only if e*^ = ib|c|d/c|d|. Plugging these values into 
h gives the points where the length ol g{9) is changing most and least with 
respect to 9. 

To find the radius of T', one computes /(27r)/27r: 

1{2tt) 1 /"^^ \ad-bc\ 

-da 



27r 27r Jq |c|2 + |(i|2 + cde^e + ^de- 
\ad — bc\ f dz 



2m J {c + dz){d + cz) 

where the contour on the last integral is the unit circle. A routine evaluation 
by the residue theorem along with (j3.1ip gives the result. The center can 
easily be computed as it is the circumcenter of the triangle formed by any 
three points on the circle, for example, zi = h{l), Z2 = h{—l), and z^ = h{i). 
The well-known formula for the circumcenter of three non-collinear points 
in the complex plane 

\zi\'^iz2 - Z3) + \z2\'^iz3 - Zi) + 12:3^(^1 - Z2) 



Zlizs - Z2) + Z2{zi - Zs) + Z3{Z2 - Zi) 

thus gives the center of the sequential closure circle. The final conclusions of 
the theorem follow immediately from the formulas for the points of highest 
and lowest concentration. □ 

Corollary 6. If the sequential closure of the continued fraction in ()3.6p is 
a line in C, then the point of highest concentration of approximants in the 
sequential closure is exactly 

h(oo) + h(0) I /a b 
2 =2[-c^-d 

the average of the first two modifications of (|3.6|) given in Corollaryl^ More- 
over, if the sequential closure is M, then the limiting probability density func- 
tion for the approximants is given by the Cauchy density function 

(3.43) p{x) = —r- TFT^' 
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where 6 is the dispersion (scale) parameter given by 
^ /i(oo) - h{0) I fa b 



2i 2i \c d 

For real argument period 1 continued fractions (not limit periodic) the 
density function has been studied before (somewhat more informally than 
here); see [331134]. 

Proof. If the sequential closure is a line, then Theorem [8] implies that |c| = 
\d\. The same theorem also implies that the point of highest concentration 
is given by 

\c\ + \d\ ■ 

When |c| = \d\, this simplifies to 

I fa b 
2\c^d 

the average of h{oo) and /i(0). 

Suppose the sequential closure is M. Let the point x £ M be related to 
the point z on the unit circle via 

, , , az + b 

X = h{z) = — , 

cz + d 

and suppose z = e*^. Let £ (0, 27r] be the angle for which z is mapped 
to oo by h{z), and put zq = e^^°. Let p{x) denote the probability density 
function and let fi denote the i-th. approximant of (j3.6p . Then for any 
interval [a,b], 

\ X , 1- #{/» £ [a,b]}o<i<n 
p[x)dx = lim 

_ l,ih-^{[a,b])) 

where, the second equality follows from remarks made in the discussion 
preceding Theorem [HI In particular, 

, , , length of the arc clockwise from zq to z Oq — 
P^'^'' = 2^ = 

Using the Fundamental Theorem of Calculus, one obtains 

, , —IdO —1 dz ad — be 

p{x) 



2tt dx 2TTiz dx 27ri{cx — a){dx — b) 
h{oo) - h{0) 6 



2m{x - h{oo)){x - /i(0)) t:[{x - xqY + 5"^] ' 
where the parameters are as given in the corollary. □ 
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Note that our arguments above (in the discussion preceeding Theorem 
[HI and the proofs of Theorem [8] and its corollary) apply to any uniformly 
distributed sequence on T and thus the first part of the following corollary 
follows. 

Corollary 7. Let {sn} be a sequence that is uniformly distributed on T and 
let h[z) = {az + h)/{cz + d) he a linear fractional transformation that maps 
T to M. Then the sequence has a Cauchy distribution on M with 

parameters xq = {a/c + b/d)/2 and 6 = {a/c — b/d)/{2i). 

Conversely, every Cauchy distribution on M arises as such a Mobius trans- 
formation of a uniformly distributed sequence on T. 

Proof. Because of the above discussions we need only to prove the converse 
direction. Let a sequence {sn} have a Cauchy distribution on M with pa- 
rameters xq and S. Then the sequence {sn/S — xo/5} has a standard Cauchy 
distribution centered at with dispersion parameter 1 on M. Then the trans- 
formation {iz + l)/{—iz + 1) maps this sequence to a uniformly distributed 
sequence on T, see |30) . Finally note that the composition of these two maps 
of the sequence {s^} is a Mobius transformation. □ 

The circular Cauchy distributions of [30] are special cases of the family 
of distribution functions on circles in C obtained by applying an arbitrary 
non-trivial Mobius transformation to a uniform distribution on the unit cir- 
cle. We call this family the Mobius-Cauchy distribution family. It is clearly 
closed under the full Mobius group, unlike the circular Cauchy distributions 
which are closed under only those Mobius maps that fix the unit circle, 
[30] . As (real) Cauchy distributions are closed under the real Mobius group, 
the Mobius-Cauchy family forms a natural generalization of the usual real 
Cauchy distribution that also includes the family of circular Cauchy distri- 
butions special case. 

3.2. Computing subsequences of approximants converging to any 
point on the sequential closure. We recall one of the main conclusions 
of Theorem [71 Namely, that if |pn| < oo, X^k^l < l"^! = 1/^1 = 1 
and A = a//3 is not a root of unity, then /„, the n-th approximant of 
K(—af5 + qi\a + (3 + Pi), satisfies 

for some o, b, c and d G C. Thus the approximants densely approach a circle 
in the complex plane and a natural question is the following: is it possible to 
explicitly determine a subsequence of approximants converging to /i(e^'^*^), 
for any 6 G [0, 1)? The following algorithm solves this problem assuming 
knowledge of the regular continued fraction for 9. 

Let A = e^'^^T, 7 £ (0,1) and let {a„/6„} denote the sequence of even 
indexed approximants in the regular continued fraction expansion of 7. Since 
A is not a root of unity, it follows that 7 is irrational. For real z, let {z) 
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denote the fractional part of z. Thus {z) = z — \_z\ . Let 9 £ [0, 1) and, for 
n > 1, let Vn denote the least positive integer satisfying < rn/bn — d < l/bn- 
For any positive integer x, 

X7 - = X f 7 - V ^^^^ + f ?^ 

V bnj bn \bn 

Since gcd(a„, 6„) = 1, there exists a non-negative integer x < bn satisfying 
a„x = r„ (mod 6„). Let kn be this solution. Since (a„/cn — fn)/bn G it 
follows that 

If the sequence {kn} is unbounded, let {j„} be a strictly increasing sub- 
sequence. If {kn} is bounded, replace each kn by kn + bn and once again 
let {in} be a strictly increasing subsequence. From the theory of regular 
continued fractions we have that, in either case. 



h 

and thus that 



On 

b 



n 



< (kn + bn) 



an 

bn 



2 

bn 



Unl -9)^0. 

It now follows that x /i(7^") h{e'^'''^). Thus 

lim = /^(e^-^). 

n— >oo 

Note that for rational A = m/n, one takes approximants in arithmetic 
progressions modulo n to obtain the subsequences tending to the discrete 
sequential closure. 

3.3. Related work. We are aware of four other places where work re- 
lated to the results of this section was given previously. Two of these 
were motivated by the identity (jl.SI) of Ramanujan. The first paper is [2] 
which gave the first proof of (jl.Sp . The proof in [2\ is particular to the 
continued fraction (jl.Sh . However, section 3 of [2j studied the recurrence 
Xn = + an-i)xn-i — Xn-2 and showed that when X^„>x |a„| < oo, the 
sequence has six limit points and that moreover a continued fraction 
whose convergents satisfies this recurrence under the h assumption tends to 
three limit points (Theorem 3.3 of [2]). The paper does not consider other 
numbers of limits, however. Moreover, the role of the sixth roots of unity 
in the recurrence is not revealed. In the section 6 of the present paper, we 
treat the general case in which recurrences can have a finite or uncountable 
number of limits. Previously in [10] we treated such recurrences with a finite 
number m > 1 of limits as well as the associated continued fractions. 

Ismail and Stanton [18] gave a proof of (jl.Sp and also obtained Corollary 
[9] below. Their method was to use properties of orthogonal polynomials 
and they obtained theorems on continued fractions with any finite num- 
ber of limits. Unfortunately the method has limitations on the perturbing 



36 



DOUGLAS BOWMAN AND JAMES MC LAUGHLIN 



sequences, requiring them to be real. On the other hand, when the ap- 
proach apphes, it yields a formula for the limits of the continued fraction 
in terms of associated orthogonality measures. The general theorem of {TS\, 
Theorem 5.2, is actually a simple application of Theorem 40, of Nevai [32j . 
The overall approach of [H] was actually employed previously in [T]. [18] 
also contains a number of other beautiful explicit new continued fraction 
evaluations, similar to (jl.Sp . 

We also mention here, in as much as it deals with the convergence of 
subsequences of approximants of continued fractions, the results obtained 
when the approximants of a continued fractions form normal families. In 
these cases there are theorems for expressing the limiting function (of a con- 
vergent subsequence of approximants) as Stieltjes integrals. See for example 
Henrici [16] and Wall [50]. The work of the present paper deals, however, 
with the pointwise limits of the continued fractions, rather than the limits 
of subsequences of functions of the variable z in a continued fraction. 

Finally, we briefly compare our results with a theorem of Scott and Wall 
[m [50] . Consider the continued fraction 

Theorem 9 (Scott and Wall). If the series ^ |^2p+i| cind J2 |^2p+iSp|, where 
Sp = b2 + bi + • ■ ■ + b2p, converge, and liminf \sp\ < oo, then the continued 
fraction ()3.44p diverges. The sequence of its odd numerator and denomina- 
tors convergents, {A2p+i} o.nd {i?2p+i}; converge to finite limits Fi and Gi, 
respectively. Moreover, if s is a finite limit point of the sequence {sp}, and 
lim Sp = s as p tends to oo over a certain sequence P of indices, then A2p 
and B2p converge to finite limits F{s) and G{s), respectively as p tends to 
oo over P, and 

FiG{s) - GiF{s) = 1. 

If the sequence {sp} has two different finite limit points s and t, then 

F{s)G{t) - F{t)G{s) = t-s. 

Finally, corresponding to values of p for which limsp = cxd, we have 

r A2p Fi 
lim — - = —- , 

B2p Gi 

finite or infinite. 

One similarity of this theorem to the present work is that it makes no 
assumptions about the size of the sequential closure. It retains much of 
the structure of the Stern-Stolz theorem, in as much as it focuses on the 
parity of the index of the approximants. However, to understand sequential 
closures in general, all subsequences need to be considered. At any rate. 
Theorem [9| does not focus on the sequential closure, but rather on loosening 
the li assumption to the subsequence odd indexed elements of the continued 
fraction. 
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4. A GENERALIZATION OF CERTAIN RAMANUJAN CONTINUED FRACTIONS 

In this section we study the non-trivial case of Theorem \7\ in which the 
perturbing sequences pn and g„ are geometric progressions tending to 0. The 
inspiration for this is the beautiful continued fraction ()1.4p of Ramanujan. 
Theorem 1101 below covers the loxodromic (convergent), parabolic, (conver- 
gent in this case), as well as the elliptic (divergent) cases simultaneously. 
Another point of this section is that it shows how Theorem [7] gives another 
approach to the problem of evaluating continued fractions. In fact it is in- 
teresting to compare the proof of Theorem [TO] to the proofs of special cases 
given previously by different methods, see [2| [TOl [T8] . 

We first recall that a basic hyper geometric series is defined for \q\ < 1 

by 

Ma; b; q,x)=f2 , (-l)"g"("-^)/^x-. 

^ {q;q)n{b; g)„ 

For the g-product notation used here, please see the introduction. Recall 
the notation defined before Theorem [71 Now define T^, the parabolic 

unitary characteristic, to be the map from C to 2^, equal to the set T^, 
when A 7^ 1, and T^^ is a set consisting of any fixed element of C — {1}. (The 
element of C — {1} does not matter.) Again, as in Theorem [71 we assume 
T*^^ annihilates inconvenient terms when \a\ ^ \f3\ and log^(a//3) G Z. Note 
that in these cases, the continued fraction is convergent and thus asymptotic 
to its limit, which will be given by the quotient of the remaining terms on 
the right-hand side. Thus in ()4.2p the restriction log^(a//3) ^ Z is not too 
serious. 

Theorem 10. Let q,a, [3 ^ C, a/(3 £C and \q\ < 1, then 

—aj3 + xq —a[3 + xq^ —a[5 + xq^ 
a + P + yq + a + p + yq^ + a + P + yq^+-- - ~ 

Moreover, assuming additionally that a ^ [i and logg(a//3) ^ Z, 

-aj3 + xq —a(5 + xq^ —a/3 + xq^ —a(3 + xq^ 



(4.2) 

a + P + yq + a + /3 + yq^ + a + (3 + yq^ + • • • + a + fi + yq 



This theorem contains an evaluation for a (/-continued fraction under the 
widest possible conditions. Notice that the inherent symmetry between a 
and fj is explicit on both sides of the equations. We present Theorem [lOl as 
a model of the evaluation of (/-continued fractions (for complex \q\ < 1). 
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In this theorem, we have not provided the error term for the difference 
between the left and right hand sides of (j4.2[) . But Theorem [7] implies 
that in the elliptic case (when \a\ = = 1), this error is 0{q^). In the 
loxodromic case (|a| 7^ |/3|), the error term can be computed from Corollary 
11 in Chapter IV of [28]. 

Comparing this theorem to Theorem [7l it is natural to enquire about 
the values of the parameters a, 6, c, and d. In fact, the proof of Theorem 
[To] follows the structure of Theorem [7] and the constants are the expressed 
functions in the above statement. 

Remark. The need for using the notion of sequential closure instead of 
the set of accumulation points can easily be seen here if one considers the 
case X = y = q = and a = —(3 = 1. The approximants of the continued 
fraction on the left hand side of (j4.ip form the sequence {00, 0, 00, 0, . . . }. 
Accumulation points are defined for sets and as a set it has no accumulation 
points. On the other hand, one cannot use the closure of a sequence, since 
in general cases (where q ^ 0) the approximants are not in the sequential 
closure. Note that in this trivial case, a/ (3 = —1 so that T*^^ = {1,-1}, 
and the right hand side of (gT]) simplifies to (T*/^ + 1)/(T*/^ - 1) so that 
both sides are equal as sets. 

Before proceeding with the proof, we note a few simple corollaries. The- 
orem [10] generalizes certain well-known continued fraction evaluations. For 
example, setting a = y = and (3 = 1, dividing by x, changing x to x/q, and 
taking reciprocals in ()4.ip yields the evaluation of the important generalized 
Rogers- Ramanuj an continued fraction: 

Corollary 8. For x,q & C and \q\ < 1, 



^ ^ xq'- _ Z^m>0 (q) 



1 + 1 + 



Z^m>0 [q)r. 



The next corollary generalizes Ramanujan's continued fraction (|1.4p with 
three limits given in the introduction. 

Corollary 9. Let uj be a primitive m-th root of unity and let u) = l/uo. Let 
1 < i < m. Then 

, . , 1 1 1 

(4.3 lim ^ J— 

_ 101 (O; gcj^; g, -g^cj) - i0i (O; qju?; q, -q^/uj) 
UJ-'' (0; g, -qu) - icpi (0; q/uj'^;q, -q/uj) 

Proof. This is immediate from (j4.2p . upon setting x = 0, y = 1, a = uj, 

(3 = uj-"^, n = mk + i, then noting that uj"^^ = 1. □ 



This result in its present form first appeared in ^18j. The authors found 
it independently and gave a different proof in [10] . 
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Corollary 10. 

2 3/2 - 3/2 - 3/2 

In fact, the limiting probability density function of the approximants of this 
continued fraction is given by 

V7 



^^^^ 27r(2x2 - 3x + 2) • 

Proof In TheoremUnitake a = 3/4+i\/7/4, /? = 3/4-i\/7/4, q = x = y 
The limiting probability density function follows from Corollary [6l 



0. 

□ 




Figure 1 . The distribution of the first 3000 approximants of 
3/2 + K'^^i with the point x = 3/4 of predicted highest 
concentration and the limiting probability density function 
p{x) = ^/7/(27r(2x2 - 3x + 2)). 

Thus Theorem [TU] unifies the celebrated Rogers- Ramanuj an continued 
fraction, Ramanujan's continued fraction with three limits, and (II. 2p . the 
continued fraction for R given in the introduction; it gives identities that 
can have any rank in Z"*" U {c}. 

It is interesting to consider that although identities such as that in Corol- 
lary [10] may appear useless as they display the oscillatory divergence of the 
continued fraction, the divergence is actually quite well controlled in the 
sense that there is an explicit probability density function for the approx- 
imants. Later Corollary [6] will be used to compute the point of highest 
concentration as the average of two convergent variants of the original con- 
tinued fraction. Thus these continued fractions come with a kind of built-in 
summability. Of course, if one is interested in computing the sequential 
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closure circle, the continued fraction converge to it rather rapidly, as was 
remarked in the first paragraph following Theorem [101 See also Figure [2] 
below. 

Here is the proof of Theorem [101 Following the proof other special cases 
are studied. 

Proof. This theorem is a simple consequence of Theorem [7] and work from 
our paper pjj. First consider the case \a\ = = 1. Let Pn and Qn denote 
the nth numerator and denominator convergents of (j4.2jl . In [llj we worked 
with the related continued fraction 

1 —a(5 + xq —afi + xq^ —af) + xq^ 

l + a + l3 + yq + a + (3 + yq^ + a + (3 + yq^^ ' 

(Note that in [11] a, 6, c, and d were written for a, /3, x, and y, respectively.) 
Let An and Bn denote the n-th numerator convergent and n-th denominator 
convergent, respectively, of (|4.4I) . Then it is immediate that Pn = Bn+i — 
An+i and Qn = An+i- Observe that since both of these continued fractions 
are symmetric in a and /?, the explicit expressions for An and Bn ((2.5-6) 
from [TT]) are also valid with a and (3 interchanged. All that is needed is to 
compute the limits (j3.8p from Theorem [71 Now employing the convergents 
of (j4.2p instead of those for ()4.4p puts equations (2.15) and (2.16) from [11] 
into the forms 



i2\ " 



(AK^ Pn _ (xq/P-a) ^ {-xq/yP)n . n(n.-i)/2 ( -y<l 

^^•^^ A^™oo/3^ l-a//3 A. (ag//3)„(<;>/ « \ (i } 

and 

f4 6) lim i V {-Xqly(i)n .n n(n-X)l2 ( 'V^Y 

^^•^^ N^ooH^ l-a/(3^^{aq/P)n{q)n^- ' ^ \ P ) 

Thus, for example, using (13. Sh . 



h = - lim I3-^{Pn - aPN 



{ lim (3-^Pn - (a/p) lim /J-^^^^^Ptv-i 



(1-Q//3) lim 13-^Pn 



N 

2 



xq \ , I ^xq aq —yq 

a 101 



(3 J '^'\yf3' (5 (3 

Here the last three equalities followed from (j4.5p . 

The identification of the other constants in h is similar, except that one 
must interchange the role of a and (3 in (14. 5p and (j4.6p when calculating a 
and c. 

The case |a| = \f3\ ^ follows by taking the equivalence transformation 
as in the proof of Theorem [71 Note that the parameters a and (3 in the 
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4> function are not restricted to this domain. For brevity in the rest of 
the proof, we employ the variables a, 6, c, and d, with the understanding 
that they represent the cj) functions with the above arguments, analytically 
continued to their widest domain. (It is easy to check that the (j) function is 
meromorphic in its variables in the complex plane.) 

Now assume that \a\ ^ \(3\. First note that the difference equation 

(4.7) y„ = (1 + A - V)^n+i + (-A + s V)y„+2 

has a solution y„ = i0i(s; Ag; g, zg"). (This can be checked simply by 
equating coefficients.) By Auric's theorem, see Corollary 11, Chapter IV of 
[28], this solution of (j4.7[) is minimal if | A| < 1 or A = 1, and thus for |A| < 1 
or A = 1, 

i(t)i{s\ Xq; q, z) _^ , , , -\ + sz -\ + szq 

L + A — Z + 



1^1 (s; Xq; q, zq) 1 + X — zq + 1 + A — zq"^ + • • • 

Putting s = —(3~^xy^^q, X = a/(3, and z = —P~^yq, taking reciprocals, 
multiplying both sides by —a + xq/ (3 and applying a simple equivalence 
transformation to the continued fraction, yields that for \a\ < \(3\ or a = 

b —a(3 + xq —aj3 + xq^ —a(5 + xq^ 
d~a + P + yq + a + /3 + yq'^ + a + P + yq^-\ ' 

For \a\ > \P\, symmetry gives that 

a —a(5 + xq —af5 + xq^ —a(3 + xq^ 
c~a + p + yq + a + p + yq^ + a + f3 + yq^+---' 

The conclusion follows by noting that for \a\ < |/?| as n — > oo, 

aA"+^ + b _ b 
cA"+i + d ^ d' 

while for |q| > 

aA"+^ + 6 a 



cA"+i +d c ■ 

□ 

Remark. We could have simply used (|4.5|) and (|4.6|) to complete the 
\a\ / \(3\ part of the proof, but the approach via Auric's theorem seemed 
preferable as it also yields the evaluation of the continued fraction in the 
parabolic case a = /5 / 0. 

Consider the special case of the continued fraction in the theorem in which 
X = and y = 1- Then 

-P i<^i (0;^;q,^)z + a i<Ai (o; f ; g, ^ 
(4.10) h{z) = ^ ^ ^ 



i'?'UO;§;'z,^U- i<Ai (o;f;g,f 
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and thus that the sequential closure of the continued fraction 

-a/3 a(3 af3 
G[a,p,q) := ^ 5- 

is on the circle /(T), where / is defined by 




Figure 2. The convergence of G(exp(zA/7), exp(zA/5), 0.1) 



Figure E] shows the first 3100 approximants of G(exp(i\/7), exp(^\/5), 0.1) 
and the corresponding circle f{T) predicted by the theory. The larger dots 
show the points, again predicted by the theory, of highest and lowest concen- 
tration of approximants. Note the error, e„ = O(10~") and experimentally, 
min^gf \^n/Bn — f{z)\ ~ 10~" in agreement with the theory. 

Figure [3] shows the first 2700 approximants (with the first 55 omitted, 
since they lay relatively distant from the circle of convergence) of 
G'(exp(z^/7),exp(z(A/7 + 27r/ll)),0.99exp(iA/T7)) and its convergence to the 
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eleven limit points /(2/c7r/ll), < A: < 10, where f{z) is the associated 
linear fractional transformation, together with the circle /(T). The value 
q = 0.99exp(^\/T7)) was chosen to be close to 1 in absolute value, with 
the aim of slowing down the convergence so as to make the behavior more 
visible. 
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The Bauer-Muir transformation (Proposition [2]) , can applied to Theorem 
[To] to obtain convergent continued fractions. The results are contained in 
the following corollary. The continued fraction from Theorem [10] will be 
denoted as follows. 

—aP + xq —aj3 + xq^ —aj3 + xq^ 

K{a, 13, X, y) := — — — — — — ^ — — — 3 

a + j3 + yq + a + p + yq'^ + a + fJ + yq-^ + ••• 

Corollary 11. Let \q\ < 1 and let \a\ = |/3| 7^ be distinct complex numbers 
such that arg(a//3) is not a rational multiple o/vr. When the approximants 
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of K(a, (3, X, y) are real, so that the continued fraction is dense on M, its 
point of highest concentration is given by 



xo 



If f3yq + xq ayq + xq 

-pH —, — r-a + 



2 V a + yq + K{a,l3q,xq,yq) f3 + yq + K{aq, P,xq,yq) 

In fact, the approximants of K{a, (3, x, y) have the Cauchy distribution 

6 



p{x) 



7r[(x — xqY + ^'^] ' 



with scale factor 

1 / Pyq + xq ayq + xq 

= — { -p -\ ■ 777 — S 7 + « 



2i V a + yq + K{a,f3q,xq,yq) f3 + yq + K{aq, (3,xq,yq) 

Moreover, assuming only the first sentence, 
(4.11) 

f xq R\ A\ ( ~y'P 1 

\— - 1^) Pyq + xq 

h[00) = -: ^ r = —p + 



and, 
(4.12) 



^(p^y-^-§R-q^^] a + yq + K{a,Pq,xq,yq) 



-0- - a] icpi I -0-; q, —0- l ^yq _)_ xq 



^^^^ ^(^^ "''^ P + yq + K{aq,(3,xq,yq)' 

This corollary exhibits the interesting phenomenon of a continued fraction 
that diverges (and is dense in R), yet the statistical distribution is given as 
a rational function of the same continued fraction evaluated at points where 
it converges. 

Proof. By Corollary [3l it follows that h{oo) = lim„_^oo /n(— /3) and h{0) = 
lim^^oo /n(— «)) where fn{w) denotes the nth modified approximant of 
K{a, (3,x,y). Apply Proposition [2] to the continued fraction K{a, (3,x,y) 
with the modifications Wn taken to be the constant sequences {— /?, — /3, . . . } 
and {—a,—a,...}, respectively, to obtain the convergent continued frac- 
tions for /i(oo) and /i(0). The series expressions for /i(oo) and /i(0) follow 
directly from Theorem [TOj The final conclusion is then an application of 
Corollary [6l The conditions x ^ —fiy and x ^ —ay are needed to ensure 
that the Bauer-Muir transformations exist (A„ 7^ in Proposition [2]). Ob- 
serve that these conditions may be dropped in the corollary by the following 
well-known version of the g-binomial theorem [15] : 



n=0 



□ 
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Notice that the two continued fraction to series identities in this corollary 
are equivalent by the symmetry in a and /?. Also, the series equals continued 
fraction identities follow from Theorem [10] by setting /? to and then 
assuming \a\ = 7^ 0, and manipulation of the continued fraction and 
series. (A contiguous relation needs to be employed to simplify the series in 
this case.) 

In some cases the infinite series in the corollary above can be expressed 
as infinite products. 

Corollary 12. Let \q\ < 1. Then 

{q;q^)oo{q^;q^) 



00 



Proof. In (j4.1ip . set x = 0, y = 1, replace q by g^, set (3 = —q and a = 
q and simplify the resulting continued fraction by applying a sequence of 
equivalence transformations. 

For the right side we use two identities due to Rogers [30] (see also 
and IS], identities A. 16 and A. 20): 



n(n+2) 



„=0 (9^ 9^)00(93; g5)^(_g2.g2)^ 



^^olg^;*?^)™ iq;q')ooiq^;q')ooi-q^;q^: 



00 



Finally, cancel a factor of q on each side □ 

The continued fraction above is clearly a transformed version of the 
Rogers- Ramanuj an continued fraction since they converge to the same limit. 

5. Applications to (r, s)-matrix continued fractions 

Levrie and Bultheel defined a generalization of continued fractions called 
(r, s)-matrix continued fractions [27]. This generalization unifies a number 
of generalizations of continued fractions including "generalized (vector val- 
ued) continued fractions" and continued fractions obtained by composing 
together higher dimensional Mobius transformations, see [25] and [S] for 
definitions. 

Here we show that our results apply to limit periodic (r, s)-matrix con- 
tinued fractions with eigenvalues of equal magnitude yielding estimates for 
the asymptotics of their approximants. 

For consistency we closely follow the notation used in [27] to define (r, s)- 
matrix continued fractions. Let OJts,r(C) denote the set of s x r matrices 
over the complex numbers. Let 0^ be a sequence of n x n matrices over 
C. Assume that r + s = n. A (r, s)-matrix continued fraction is associated 
with a recurrence system of the form = Y^-i^fc- The continued fraction 
is defined by its sequence of approximants. These are sequences of s x r 
matrices defined in the following manner. 
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Define the function / : D G 9Jt„(C) ^ ^s,r{C) by 

(5.1) f{D) = B-'A, 

wliere B is the s x s submatrix consisting of the last s elements from both 
the rows and columns of D, and A is the s x r submatrix consisting of the 
first r elements from the last s rows of D. 

Then the k-th approximant of the (r, s)-matrix continued fraction associ- 
ated with the sequence 6^ is defined to be 

(5.2) Sfc :=/(efc^fc-i---Mi). 

To apply Theorem[5]to this situation, we endow 9Jtsxr(C) with a metric by 
letting the distance function for two such matrices be the maximum absolute 
value of the respective differences of corresponding pairs of elements. Then 
when / is continuous, our theorem can be applied. (Note that / will be 
continuous provided that it exists.) 

Let lim^^oo = 6, for some 6 G 9Jl„(C). Then the recurrence system is 
said to be of Poincare type and the (r, s)-matrix continued fraction is called 
limit periodic. 

After this definition Theorem H] can be applied and the following theorem 
results. 

Theorem 11. Suppose that the condition X]^>]^ ^|| < oo holds, that the 
matrix 9 is diagonalizable, and that the eigenvalues of 6 are all of magnitude 
1. Then the kth approximant s/. has the asymptotic formula 

(5.3) Sk X fie'^F), 

where F is the matrix defined by the convergent product 

F := lim Q-^OkOk-i- ■ ■ 0261. 

fc— >oo 

Note that because of the way that (r, s)-matrix continued fractions are 
defined, we have taken products in the reverse order than the rest of the 
paper. 

As a consequence of this asymptotic, the sequential closure can be deter- 
mined from 

©(^fc) = ©{f{O^F))- 
In one general case, detailed in the following theorem, we actually get a 
convergence theorem. 

Theorem 12. Let 9k he a sequence of n x n matrices over C satisfying 

J2\\0k - 0\\ < (x^, 

k>l 

where 9 is a diagonal (or antidiagonal) matrix with all diagonal (or antidi- 
agonal) elements of absolute value 1. Let r and s he positive integers with 
r + s = n. 
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Then the matrix 



F := lim e^'^OkOk-i ■ ■ ■ 0261 

k^oo 



exists. Suppose further that the bottom right s x s submatrix of F is non- 
singular. Then the {r, s) -matrix continued fraction defined by equation (|5.2p 
converges to f{F). If 6 is antidiagonal, then the even approximants of the 
{r, s)-matrix continued fraction defined by equation ()5.2p tend to f{F), while 
its odd approximants tend to f{AF), where A is the antidiagonal matrix with 
Is along its antidiagonal. 

Proof. The matrix F exists by Theorem [J] (or more precisely, the "trans- 
posed" version of Theorem [5]) . Let 



6 = diag(Ai 



I An)- 



By (I53D, 

sk ^ fie'F) 



n— s+1, n— s+1 



F 



-1 



n— s+l,n 



... Ay V ^n,n-s+l 

\ / Fn^s^i^i 



^n-s+1 • 







AJ , 



Fn 

F 



\ Fn,l 



n— s+l,r 



7? 



^ -^n— s+1, s+1 • • • -^n— s+l,n\ / -^a— s+1,1 • • • F^ 



-s+l,r 



n,n— s+1 



n,l 



F 



n,r 



f{F). 



Thus Sk converges to the final matrix product above. 

For the case where 9 is an antidiagonal matrix, ^■^'^ is a diagonal matrix 
and the proof for the even approximants is virtually the same as for the case 
where is a diagonal matrix. If is an antidiagonal matrix, ^^fc+i -g ^^^^ 
an antidiagonal matrix. Once again by ()5.3p . 

S2k+i ^ fiO^'^'F) 
// 



^^^^"'^)ri~s+l,s\ / -^l,n-s+l • • • -^1, ^ ^ 



V Fs^n- 



-s+1 



/ ... (02fc+l) 



n— s+1, 



/ Fi^i . . . Fi^r 



32fc+l^ 



n,l 



V Fs,i 
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= fiAF), 

where A is the antidiagonal matrix with I's along the antidiagonal. Thus 
S2k+i converges to the final matrix product above. □ 

Consider now the n = 2 antidiagonal case of Theorem [T2l The matrix 6 
then has the form 



1 

1 



Choose 6k to have the form 



bk V 

Using the correspondence between matrices and continued fractions ()3.ip . 
we at once obtain the following corollary, first given in [lOj . 



Ok 



Corollary 13. Let the sequences {a„} and {bn} satisfy an ^ —1 for n > 1, 

^ |o„| < oo and '^\bn\ < oo. Then 

n=l bn 



Urn Q2n+p = -Bp / oo, 



diverges. In fact, for p = 0, 1, 
lim P2n+p = Ap^oo, 
and 

oo 

AiBo-A^Bi = J](l + a„). 

n=l 

In fact, Corollary [Ql is also the a = 1, /3 = — 1 (so m = 2), g„ = a„ and 
Pn = bn case of Corollary [S] When a„ = 0, this corollary reduces to the 
famous Stern-Stolz theorem discussed in the introduction. 

One of the main results of the paper [10] was Corollary O which we 
applied to obtain an infinite sequence of theorems, similar to the Stern- 
Stolz theorem, but with continued fractions of different ranks. Notice that 
Theorem [12] provides yet another family of generalizations. 

It is interesting to compare Corollary [13] with the "The General Stern- 
Stolz Theorem" from [5j in the case of continued fractions. The corollary 
for the case of complex continued fractions is: 
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Corollary 14. [Corollary 7.5 of \^] If |1 — |a„|| and Y2n converge, 
then K{an\bn) is strongly divergent. 

The first condition in this result is weaker than analogous condition in 
Corollary [13] above. But it should be remarked that Theorem [H Corollary 
[T3| and Corollary 1141 are, in fact, equivalent; the two corollaries follow from 
Theorem [T] by an equivalence transformation (and a little analysis). Next, 
the condition on the partial numerators in Corollary 1 141 encodes the informa- 
tion that the matrices representing the continued fraction are a perturbation 
of unitary matrices. We could have obtained the same result by using Theo- 
rem [3l however in this situation one does not obtain as detailed information 
about the limits of the convergents. In particular. Corollary 1131 also proves 
the convergence of the subsequences of convergents {Pn} and {Qn} of equal 
parity. Corollary [T3] does not furnish this part of the conclusion. On the 
other hand, it does prove strong divergence, defined in section 2. Indeed, 
the continued fraction in Corollarvll4lis not necessarily limit periodic. 



6. POINCARE TYPE RECURRENCE RELATIONS WITH CHARACTERISTIC 
ROOTS ON THE UNIT CIRCLE 

Let the complex sequence {xn}n>o have the initial values xq, . . . , Xp_i 
and subsequently be defined by 

p-i 

(6.1) X}2+p — ^ ^ Q"n,rXn+r i 

r=0 

for n > 0. It is assumed that for n sufficiently large a„^oOn,p-i 7^ 0. Suppose 
also that there are numbers oq, . . . , Op-i such that 

(6.2) lim anr = dr, < r < p — 1. 

n—foo ' 

A recurrence of the form (16. Ih satisfying the condition (16. 2p is called a 
Poincare-type recurrence. Such recurrences were initially studied by 
Poincare and later Perron who proved the Poincare- Perron theorem [24^ [37] : 

// the roots of the characteristic equation 

(6.3) tP - ap^itP-^ - ap_2t^"^ ao = 

have distinct TioTTfiSj then linij2 >oq Xn+i/xn = a, where a is a root of ()6.3p . 

Moreover, for each root a of (16. 3p . there exists a solution of (j6.ip with term 
ratio tending to a. 

Because the roots are also the eigenvalues of the associated companion 
matrix, they are also referred to as the eigenvalues of (16. ip . This result 
was improved by O. Perron, who obtained a number of theorems about 
the limiting asymptotics of such recurrence sequences. Perron [35] made 
a significant advance in 1921 when he proved the following theorem which 
treated cases of eigenvalues which repeat or are of equal norm. 
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Proposition 3. Let the sequence {xn}n>o be defined by initial values xq, 
. . . , Xp-i and by Ii6.1\) for n > 0. Suppose also that there are numbers 
gq, . . . , Cp-i satisfying ( fg.^j) . Let qi, q2, ■ ■ ■ Qa be the distinct moduli of the 
roots of the characteristic equation lid. 3\) and let l\ be the number of roots 
whose modulus is q\, multiple roots counted according to multiplicity, so that 



Then, provided Onfi be different from zero for n > 0, the difference equation 
()6.ip has a fundamental system of solutions, which fall into a classes, such 
that, for the solutions of the X-th class and their linear combinations. 



The number of solutions of the X-th class is Ix. 

Thus when ah of the characteristic roots have norm 1, this theorem gives 
that 



Another related paper is [23] where the authors study products of ma- 
trices and give a sufficient condition for their boundedness. This is then 
used to study "equimodular" hmit periodic continued fractions, which are 
Umit periodic continued fractions in which the characteristic roots of the 
associated 2x2 matrices are all equal in modulus. (Thus they are exactly 
the class of limit periodic continued fractions of elliptic type.) The matrix 
theorem in [23j can also be used to obtain results about the boundedness of 
recurrence sequences. Theorem [T3l below applies to equimodular recurrences 
as well. 

More recent is the work of R.J. Kooman [24^ [25l |26] . Kooman makes 
a detailed study of the asymptotics of Poincare type recurrences as well 
as outer composition sequences of Mobius transformations. Following our 
theorem, we compare our theorem with results of Kooman. 

Our focus is on the case where the characteristic roots are of equal mod- 
ulus but distinct. Under an li perturbation we will show that all non-trivial 
solutions of such recurrences are asymptotic to a linear recurrence with con- 
stant coefficients. Our theorem is: 

Theorem 13. Let the sequence {xn}n>o be defined by initial values xq, 
. . . , Xp_i and by (j6.ip for n > 0. Suppose also that there are numbers 
ao, . . . , ttp-i such that 



I1 + I2 + . . .l„ = p. 



limsup y| 



Xn\ = q\- 



limsup a/Ix^I = 1. 



n— >oo 



00 




< r < p — 1. 



n=0 



Put 
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Suppose further that the roots of the characteristic equation 



(6.4) 



tP - a„-itP-^ 



ap-2 



tP- 



ao = 



are distinct with values ao, . . . , ap-i of equal modulus R > 0. Then there 
exist complex numbers cq, . . . , Cp_i such that 



(6.5) 



R- 



p-i 

i=0 



O (En) ■ 



Proof. Assume first that all the characteristic roots have modulus 1. Define 





( 


Op-l 


ap-2 • 


ai 








1 





. 





M := 







1 


. 







\ 








. 1 


0/ 



By the correspondence between polynomials and companion matrices, the 
eigenvalues of M are ao, . . . , Op-i, so that M is diagonalizable. For n > 1, 
define 

/ an-i,p-i an-i,p-2 • • • an-1,1 an-i,o\ 



Dr,. :-- 





1 













\ ... 1 y 

Thus the matrices M and Dn satisfy the conditions of Theorem HI From 
(I6TD it follows that 



^ ^n+p— 1^ 
•^n+p— 2 



Xp-2 

V ^0 y 



Let F have the same meaning as in Theorem [H Part (i) then gives that 



V 



Xn+p- 
Xn+p—2 



A 



/ Xp-l\ 
Xp-2 

\ ^0 y 



()6.5p follows immediately by considering the bottom entry on the left side. 
The case of modulus R follows by renormalization. This completes the 
proof. □ 



The following corollary, proved in [lOj . is immediate. 
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Corollary 15. Let the sequence {xn}n>o be defined by initial values xq, 
. . . , Xp-i as well as (16. ip for n > 0. Suppose also that there are numbers 
Oo, . . . , flp-i such that 

oo 

\ar — an,r\ < OO, < r < p — 1. 

n=0 

Assume that the roots of the characteristic equation 

tP - ap^itP-^ - ap_2tP~^ ao = 

are distinct roots of unity oq, . . . , Op-i- Let m he the least positive integer 
such that, for all j G {0, 1, . . . , p — 1}, a™ = 1. Then, for < j < m — 1, the 
subsequence {xmn+j}'^=Q converges. Set Ij = liuin-, oo Xnm+j , for integers 
j > 0. Then the (periodic) sequence {Ij} satisfies the recurrence relation 

p-i 

In+p — ^ ^ Orln+r ; 
r=0 

and thus there exist constants cq, • • • ,Cp_i such that 

p-i 

In — ^ ^ CjC^i • 

We close this section by comparing our result with those from Kooman 
[2"6j . Proposition 1.7 from [26j appears to be most closely related to Theorem 
[T3l Kooman also gives a result of Evgrafov |14J which is also similar to 
Theorem [13) 

Consider the linear recurrence (j6.ip where 

oo 

\ar — an,r\ < OO, < r < p — I. 

n=0 

If the characteristic polynomial (j6.3p has zeros uq, . . . , a^-i with < |ao| ^ 

• • • < l^p-il; then (16. ip /las solutions u^n = + o(l)). 

One difference between Evgrafov's result and TheoremllSl is that the later 
gives an error term. However, Evgrafov's theorem does not require distinct 
characteristic roots. Kooman obtained a result generalizing Evgrafov's, and 
containing an error term: 

(Proposition 1.7 of t26jj Let ao,...,ap-i be non-zero, not necessarily 
distinct numbers with \oiq\ < ■ ■ ■ < |ap-i| and let (3 : N ^ M>o be a function 
such that limn-too (3 (n) = 0, J2'^=o Pi^) < °o, and < max|aj/aj+i| < 
liminf(/3(n + l)/f3{n)) < 1 where the maximum is taken over those i such 
that \ai\ 7^ |aj+i|. Let Kn be matrices with \\Kn\\ = 0{(3{n)). The matrix 
recurrence 



(diag(ao, ■ ■ ■ , "p-i) + L[n)xn = Xn+i (n G N) 
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has solutions {a;„ } with 

x»=<eJl + o(f;/3(/i) 

V \h=n 

fori = 0,...,p-l. 

Here Cj is the ith canonical basis element for C^. The asymptotic in 
Theorem [13] can be obtained from Kooman's Proposition 1.7 as a special 
case, but with a different error term. Kooman's error term can be weaker 
because of the assumption that 

limmf \, . > 0. 

(Note that to obtain the result corresponding to the asymptotic in The- 
orem [13] from Kooman's Proposition 1.7, one must diagonalize the limit- 
ing companion matrix, M: set CMC~^ = diag(ao, . . . , Qp-i). Then put 
Kn = CDnC~^ — diag(ao, . . . ,ap_i). One easily checks that the sequence 
Kn satisfies the conditions of Kooman's proposition.) 

7. Conclusion 

We have studied convergent subsequences of approximants of complex 
continued fractions and generalizations. There is an interesting pattern of 
relationships between the limits and asymptotics of subsequences and the 
modified approximants of the original sequence. This suggests the general 
question of in which other situations do "similar" patterns of relationships 
exist? In section 2, it was shown that (at least some of) this behavior 
extends to the setting of products of invertible elements in Banach algebras. 
More generally, are there other classes of sequences that diverge by oscilla- 
tion, but for which "nice" asymptotics for the sequences exist, thus enabling 
the computation of the sequential closure? When can the probability den- 
sity functions be computed? Even more generally, when "nice" asymptotics 
do not exist, is the sequential closure interesting or useful? 

Acknowledgements. The authors would like to thank the following people 
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